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Introduction 


The  coastal  zone  is  complicated  by  internal  boundary  layers  and  growing  waves  in  offshore 
flow,  incoming  swell  and  shoaling,  and  diumally  varying  wind  fields.  The  case  of  offshore  flow 
of  warm  air  over  cooler  water  (typically  occurring  in  the  afternoon)  has  alluded  modeling 
efforts  to  date,  interpretation  of  data  analysis  is  made  difficult  by  the  fact  that  different 
influences  are  interrelated  so  that  one  can  find  significant,  misleading,  correlation  through 
indirect  relationships.  For  example,  the  roughness  length  appears  to  be  a  strong  function  of 
wave  age  in  offshore  flow.  However,  wave  age  is  strongly  correlated  with  fetch,  as  is  the 
strong  turbulence  advected  from  land,  which  decays  in  the  downstream  direction.  Based  on 
statistical  analyses,  one  could  highlight  wave  state  as  the  dominate  influence  even  though  its 
effect  is  secondary  and  happens  to  correlate  with  the  more  direct  influence  of  advection  of 
turbulence  on  the  surface  stress. 

With  this  warning,  the  present  research  has  identified  a  number  of  solid  influences.  Shallow 
convective  internal  boundary  layers  in  offshore  flow  suppress  well  developed  convective 
eddies  so  that  the  heat  flux  is  less  than  predicted  by  existing  formulations;  that  is,  the  thin 
boundary-layer  depth  influences  the  flux-gradient  relationship  at  the  surface,  requiring 
generalization  of  Monin-Obukhov  similarity  theory.  The  nature  of  the  turbulence  over  land 
exerts  a  strong  influence  on  both  the  boundary  layer  structure  over  the  sea  and  the  surface 
stress  and  fluxes  through  advection  of  turbulence  energy.  These  effects  are  significant  in  the 
near  coastal  zone,  typically  the  first  5  Or  1 0  km  offshore.  Thirdly,  confused  multimodel  seas 
without  a  dominate  spectral  peak  is  a  better  indicator  of  laYge  stress  (large  roughness  length 
and  large  Chamock  coefficient)  than  wave  age. 

On  the  other  hand,  Monin-Obukhov  similarity  theory  was  a  good  approximation  for  the  flux- 
gradient  relationship  for  onshore  flow.  In  general,  there  is  no  fundamental  difficulty  with 
Monin-Obukhov  similarity  for  stationary  homogeneous  flow  over  the  sea,  in  contrast  to  some 
earlier  suspicions. 

The  following  report  consists  of  three  papers  from  1999,  all  based  on  data  from  the  Risoe  Air 
Sea  Experiment  (RASEX).  “Observations  of  non-dimensional  wind  shear  in  the  coastal  zone” 
appeared  in  the  Quarterly  Journal  of  Royal  Meteorological  Society  (vol.  128,  pp  2685-2702). 
“Boundary-layer  transitions  in  offshore  flow  has  been  recently  submitted  to  Boundary-Layer 
Meteorology.  “The  coastal  zone”  recently  appeared  as  a  chapter  in  the  book,  Air-Sea 
Exchange:  Physics,  Chemistry  and  Dynamics,  Kluwer  Academic  Publishers,  1999. 
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19  Sept  1999 


Abstract 


The  adjustment  of  the  boundary  layer  immediately  downstream  from  a 
coastline  is  examined  in  terms  of  2  levels  of  eddy  correlation  data  collected 
on  a  mast  at  the  shore  and  6  levels  of  eddy  correlation  data  collected  from  a 
mast  2  km  offshore  during  RASEX.  The  characteristics  of  offshore  flow  are 
studied  in  terms  of  case  studies  and  inter-variable  relationships  for  the  entire 
one  month  data  set.  A  turbulent  kinetic  energy  budget  is  constructed  for 
each  case  study. 

Although  the  buoyancy-generation  of  turbulence  is  generally  small  com¬ 
pared  to  shear-generation,  the  weakly  stable  and  weakly  unstable  cases  ex¬ 
hibit  completely  different  vertical  structure.  With  flow  of  warm  air  from  land 
over  cooler  water,  modest  buoyancy  destruction  of  turbulence  and  reduced 
shear-generation  of  turbulence  over  the  less-rough  sea  surface  cause  the  tur¬ 
bulence  to  weaken  downstream  from  the  coast.  Shear-generation  at  higher 
levels  and  advection  of  stronger  turbulence  from  land  lead  to  a  maximum  of 
stress  and  turbulence  energy  above  the  surface,  in  contrast  to  conventional 
boundary  layers  where  the  stress  is  a  maximum  at  the  surface  and  decreases 
with  height. 

With  flow  of  cool  air  over  a  warmer  sea  surface,  a  convective  internal 
boundary  layer  develops  downstream  from  the  coast.  However,  an  overlying 
thick  layer  of  downward  buoyancy  flux  (virtual  temperature  flux)  is  main¬ 
tained  by  shear-generation  in  the  accelerating  offshore  flow. 
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1.  Introduction 


The  response  of  the  atmosphere  to  surface  discontinuities  is  often  posed 
in  terms  of  internal  boundary  layers  (Garratt,  1990).  In  flow  of  warm  air 
from  a  rough  land  surface  over  a  cooler  sea  surface,  the  flow  near  the  surface 
may  not  adjust  to  the  new  surface  in  the  usual  sense.  The  flow  above  the  thin 
stable  surface  layer  (Figure  la),  which  was  part  of  the  boundary  layer  over 
land,  may  become  partially  decoupled  from  the  surface,  accelerate  and  form 
a  low  level  jet  (Smedman  et  al.,  1995;  Tjemstrom  and  Smedman,  1993).  The 
shear  on  the  underside  of  the  jet  may  eventually  generate  turbulence  and 
reestablish  a  surface-based  boundary  layer.  In  this  case,  the  main  source  of 
turbulence  is  elevated  and  not  at  the  surface  and  Monin-Obukhov  scaling  is 
invalid  (Smedman  et  al.,  1995).  Further  downstream,  the  flow  can  become 
near  neutral  as  the  air  cools  and  substantially  reduces  the  air-sea  temperature 
difference  (Smedman  et  al.,  1997). 

In  some  cases  of  flow  of  warm  air  over  a  cooler  surface,  a  well-defined 
stable  internal  boundary  layer  develops  (Garratt  and  Ryan,  1989).  In  a  nu¬ 
merical  study  of  the  influence  of  the  continental  diurnal  variation  on  offshore 
flow,  Garratt  (1987)  found  that  the  onset  of  daytime  convective  turbulence 
was  advected  offshore  as  a  jump  that  could  be  traced  for  hundreds  of  kilo¬ 
meters  offshore,  well  beyond  the  fetch  in  the  present  data.  The  numerical 
simulations  of  Mengelkamp  (1991)  indicate  that  the  top  of  the  stable  internal 
boundary  layer  can  be  defined  in  terms  of  a  minimum  in  the  vertical  profile 
of  turbulence  kinetic  energy  while  numerical  simulations  of  Garratt  (1987) 
indicate  a  minimum  of  eddy  diffusivity  at  the  top  of  the  stable  internal  bound¬ 
ary  layer.  The  level  of  minimum  turbulence  separates  the  underlying  internal 
boundary  layer  from  overlying  decaying  turbulence.  In  Mengelkamp  (1991), 
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this  overlying  decaying  turbulence  still  exhibits  some  upward  buoyancy  flux 
offshore,  characteristic  of  the  upstream  convective  boundary  layer. 

A  number  of  studies  document  development  of  well-defined  convective 
internal  boundary  layers  in  flow  of  cool  air  over  a  warm  surface  (see  review 
in  Garratt,  1990).  The  growing  convective  internal  boundary  layer  is  usually 
capped  by  a  thin  entrainment  zone  of  downward  buoyancy  flux.  In  contrast, 
Sun  et  al.  (1998)  study  a  convective  internal  boundary  layer  which  is  capped 
by  a  relatively  thick  layer  of  downward  buoyancy  flux,  maintained  by  elevated 
shear-generation. 

The  above  studies  suggest  complex  and  more  varied  vertical  structure  of 
the  convective  internal  boundary  layer  in  offshore  flow.  Sun  et  al.  (2000) 
show  that  close  to  the  coast,  advection  from  land  dominates  the  near-surface 
stress.  We  will  refer  to  this  region  as  the  “near-coastal”  zone  (Figure  1). 
In  offshore  flow  over  the  near-coastal  zone,  existing  similarity  theory  may 
be  invalid  (Vickers  and  Mahrt,  1999).  In  contrast,  Monin-Obukhov  scaling 
successfully  describes  the  turbulence  energy  budget  (Edson  and  Fairall.,  1998; 
Wilczak  et  al,  1999)  and  the  flux-gradient  relationship  (Vickers  and  Mahrt, 
1999)  above  the  wave  boundary  layer  for  stationary  onshore  flow. 

Downstream  from  the  near-coastal  zone,  the  stress  is  dominated  by  wave 
growth  and  advection  of  turbulence  is  no  longer  important.  However,  the 
atmospheric  boundary-layer  depth  and  vertical  structure  may  still  be  adjust¬ 
ing.  We  refer  to  this  region  as  the  “outer  coastal  zone”.  This  partitioning 
into  two  zones  is  somewhat  arbitrary  and  over-simplified,  and  the  various 
adjustments  are  gradual  without  a  sharp  transition  in  the  downstream  di¬ 
rection.  In  fact,  definition  of  an  internal  boundary  layer  will  not  be  obvious 
from  the  observations  of  stable  offshore  flows.  However,  this  partitioning  will 
help  organize  our  initial  investigations  and  provide  some  link  to  historical 
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thinking. 

This  study  emphasizes  the  near-coastal  zone  and  examines  the  vertical 
structure  of  both  stable  and  unstable  offshore  flow  with  a  heavily  instru¬ 
mented  offshore  tower,  described  in  the  next  section. 

2.  Data 

We  analyze  offshore  tower  data  collected  during  the  Ris0  Air  Sea  Exper¬ 
iment  (RASEX).  The  full  RASEX  experiment  instrumentation  is  described 
in  Barthelmie  et  al.  (1994)  and  H0jstrup  et  al.  (1997).  In  this  study,  we 
analyze  observations  taken  at  the  sea  mast  west  tower,  located  2  km  off  the 
northwestern  coast  of  the  island  of  Lolland,  Denmark,  in  4  m  of  water,  for 
the  intensive  observing  period  3  October  through  8  November  1994.  The 
variation  in  mean  water  depth  due  to  tides  is  only  about  0.3  m.  Local  off¬ 
shore  (southerly)  flow  is  characterized  by  a  sea  fetch  ranging  between  2  km 
and  5  km.  On-shore  flow  has  a  fetch  between  15  km  and  25  km  as  it  travels 
across  an  inland  sea,  and  is  still  potentially  fetch-limited.  Fetch  is  the  dis¬ 
tance  along  the  flow  from  the  coast  to  the  sea  mast.  Water  depths  for  the 
longer  fetches  range  from  4  m  to  20  m.  The  nearby  land  surface  is  relatively 
flat. 

Various  corrections  to  the  data  are  recorded  in  Mahrt  et  al.  (1996).  Aver¬ 
aged  vertical  profiles  of  the  buoyancy  flux  and  friction  velocity  are  computed 
from  the  6  levels  of  sonic  anemometers.  In  offshore  flow,  the  vertical  vari¬ 
ation  of  the  flux  is  much  larger  than  differences  between  individual  sonic 
anemometers.  For  some  analyses,  the  fluxes  for  the  three  lowest  levels  will 
be  averaged  since  the  effect  of  instrumental  differences  may  be  larger  than 
the  actual  vertical  variation  over  such  short  vertical  distances. 
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2.1  Turbulence  energy  budget 


Insight  into  the  vertical  structure  of  the  offshore  flow  can  be  gained  by 
evaluating  the  turbulence  kinetic  energy  budget  using  the  stress  and  virtual 
temperature  flux  computed  from  sonic  anemometers  located  at  3,  6,  10,  18, 
32  and  45  m  on  the  offshore  tower. 

For  shear-generation  of  turbulence  kinetic  energy,  the  mean  wind  shear  is 
computed  from  seven  cup  anemometers  (P224b  sensor)  located  at  7,  15,  20, 
29,  38,  43  and  48  m.  Corrections  to  the  cup  anemometers  were  made  by  com¬ 
positing  the  wind  speed  profile  based  on  all  of  the  records  with  fetch  greater 
than  10  km  and  near-neutral  conditions  (abs.(z/L)  <  0.1).  This  averaged 
profile  is  fit  to  a  log-linear  height-dependence  and  percentage  corrections  for 
each  level  are  constructed  from  the  deviation  of  the  averaged  profile  from 
the  log-linear  fit.  These  corrections  partially  remove  small  systematic  irreg¬ 
ularities  in  the  profile.  Percentage  corrections  to  the  wind  speed  are  always 
less  than  2%  but  exert  a  greater  influence  on  the  shear.  The  computation  of 
the  shear-generation  term  neglects  directional  shear,  which  could  not  be  ade¬ 
quately  estimated  from  the  sonic  anemometer  data  due  to  small  uncertainties 
in  orientation. 

For  offshore  cases  (fetch  <  5 km),  advection  of  turbulence  kinetic  energy 
is  estimated  as 

yyTKEsM  —  TKElm. 

1  fetch  J 

where  V  is  the  wind  speed,  SM  refers  to  the  offshore  mast  and  LM  refers 
to  the  land  mast.  This  term  can  be  estimated  at  the  6-  and  18-m  levels, 
corresponding  to  the  common  sonic  anemometer  levels  at  the  landmast  and 
seamast.  Here,  it  is  assumed  that  the  turbulence  kinetic  energy,  TKE ,  is 
spatially  invariant  along  the  coast  for  cases  where  the  flow  was  not  perpen- 
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dicular  to  the  coast.  Unfortunately,  advection  could  not  be  estimated  above 
18  m.  Dissipation  is  estimated  following  the  spectral  approach  of  Edson  and 
Fairall  (1998).  The  spectral  slope  was  determined  from  a  least  squares  fit 
over  the  frequency  range  thought  to  be  in  the  inertial  subrange.  In  some 
cases,  the  inertial  subrange  does  not  appear  to  be  fully  developed  possibly 
due  to  nonequilibrium  conditions  in  offshore  flow  and  errors  in  the  dissipation 
estimate  may  be  significant.  We  expect  the  wave-induced  pressure  transport 
term  to  be  small  for  this  data  even  at  the  lowest  (3  m)  tower  level  (Hare  et 
al.,  1997).  The  wavelength  and  amplitude  of  the  fetch-limited  surface  waves 
axe  generally  small  compared  to  open  ocean  values. 

The  residuals  for  the  turbulence  kinetic  energy  budget  axe  expected  to 
be  laxge  because  the  pressure  transport  term  is  neglected  and  the  errors  in 
the  vertical  flux  divergence  and  horizontal  advection  of  turbulence  energy  axe 
expected  to  be  significant.  The  vertical  flux  divergence  term  (triple  correla¬ 
tion  term)  suffer  larger  random  flux  errors  compared  to  covariances,  and  the 
flux  divergence  term  is  a  small  difference  between  vertical  flux  of  turbulence 
energy  at  two  levels.  We  have  also  neglected  the  Eulerian  time-dependence 
term,  the  mean  vertical  advection  of  turbulence  kinetic  energy  as  well  as  hor¬ 
izontal  flux  divergence  of  turbulence  energy  equation.  These  terms  appear 
to  be  small  with  relatively  laxge  errors.  Considering  the  above  difficulties, 
the  residual  term  is  relatively  small  for  most  of  the  case  studies  (see  Figures 
6  and  11  below). 

2.2  Diagnostic  quantities 

To  crudely  quantify  the  unusual  vertical  structure  of  offshore  flows,  the 
stress  ratio  is  computed  as 
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stress  ratio  =  -^--upper  (2) 

stressiower 

where  stressupper  is  the  average  of  the  stress  magnitude  at  the  two  upper 
tower  levels,  32  m  and  45  m,  and  the  stressiower  is  the  average  of  the  three 
lowest  levels,  3  m,  6  m  and  10  m. 

We  compute  the  buoyancy  flux  ratio  as 

buoyancy  flux  ratio  =  hu°yancyflux»vper  /3) 

buoyancy  fluxiower 

where  the  same  levels  are  used  as  in  Eq.  2.  This  relationship  is  applied  only 
to  data  where  the  magnitude  of  the  buoyancy  flux  for  the  lower  layer  exceeds 
0.01  °Cms~1. 

2.3  External  variables 


We  now  define  two  variables  which  are  external  to  the  vertical  structure 
of  the  flow  at  the  sea  mast.  With  the  exception  of  flow  of  cold  air  over 
warmer  water,  the  turbulence  will  decay  to  some  smaller  value  downstream 
from  the  coast.  The  travel  time  of  the  flow  (fetch/wind  speed)  from  the  coast 
describes  the  period  during  which  the  decay  takes  place. 

As  a  second  external  variable,  we  define  an  “external  Richardson  number” 


D  _  hgA6/ec 

ttext  —  - n - 


0.(7 i 


(4) 


where  h  is  the  boundary-layer  depth,  crw  is  the  standard  deviation  of  the  up¬ 
stream  vertical  velocity,  A 9  is  the  difference  between  the  upstream  air  tem¬ 
perature  and  the  downstream  sea  surface  temperature,  0O  is  the  averaged 
potential  temperature,  g  is  the  acceleration  of  gravity  and  a  is  a  nondimen- 
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sional  coefficient,  arbitrarily  chosen  as  103  in  order  to  make  Rext  of  order 
one. 

With  small  values  of  Rext ,  the  strength  of  the  potential  air-surface  tem¬ 
perature  difference  is  not  sufficient  to  induce  a  well-defined  internal  boundary 
layer  in  the  face  of  vigorous  mixing  (large  crw)  due  to  the  advected  turbulence. 
In  this  case,  the  boundary  layer  advected  from  land  is  predicted  to  adjust 
to  the  new  surface  conditions  through  continuous  decrease  of  turbulence  and 
boundary-layer  depth  in  the  downstream  direction,  without  formation  of  a 
new  internal  boundary  layer  (Figure  lc).  We  refer  to  this  case  as  the  adjust¬ 
ing  boundary  layer.  While  windy  cases  of  small  air-sea  temperature  difference 
appear  to  fit  this  category,  the  RASEX  does  not  allow  examination  of  the 
spatial  variation  of  the  boundary  layer  depth.  Rext  is  external  forcing 
parameter  representing  the  influence  of  the  upstream  flow  and  sea  surface 
temperature,  as  opposed  to  the  Richardson  number  and  z  / L  over  the  sea 
which  are  functions  of  the  flow  itself  at  the  seamast. 

This  external  Richardson  number  can  be  evaluated  only  for  offshore  flow 
(fetch  =  2-5  km)  where  the  landmast  can  be  used  as  upstream  information. 
The  boundary-layer  depth  was  available  from  a  tethered  balloon  only  for 
a  subset  of  days.  We  assign  the  boundary-layer  depth  a  constant  value  of 
500  m,  in  which  case  Rext  depends  only  on  the  temperature  difference  and 
upstream  crw.  The  latter  was  estimated  as  the  average  from  the  6-  and  18-m 
levels  at  the  landmast.  The  usual  stability  parameter  z/ L  is  highly  correlated 
with  Rext  (Figure  2)  for  stable,  near  neutral  and  weakly  unstable  conditions, 
perhaps  reflecting  the  fact  that  the  local  stability  is  determined  by  advection 
from  the  upstream  land  surface.  Once  Rext  becomes  more  negative  than  -1, 
z/L  appears  to  be  independent  of  Rext,  although  there  were  no  offshore  cases 
with  large  instability. 
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3.  Land-sea  contrast 


Even  though  the  observations  are  collected  in  October  at  54  N  latitude, 
some  diurnal  signature  can  be  seen  on  many  of  the  days  with  offshore  flow. 
The  maximum  upward  buoyancy  flux  at  the  sea  surface  most  often  occurs 
in  the  morning  when  the  air  advected  from  the  land  is  coolest.  The  maxi¬ 
mum  downward  buoyancy  flux  at  the  sea  surface  tends  to  occur  in  the  late 
afternoon  when  the  air  over  land  is  warmest.  On  average,  the  diurnal  ampli¬ 
tude  of  the  buoyancy  flux  is  0.02  °Cms~1.  This  averaged  diurnal  amplitude 
is  reduced  by  the  low  sun  angle  and  numerous  cloudy  days  included  in  the 
average.  Using  temperature  at  two  levels  for  summer  for  the  same  land  and 
sea  masts,  Barthelmie  et  al.  (1996)  found,  on  average,  stable  conditions  over 
the  water  during  the  day  and  unstable  conditions  at  night. 

For  offshore  flow  for  both  stable  (R^xt  >  0)  and  unstable  conditions 
(Rext  <  0),  the  surface  stress  at  the  sea  mast  has  typically  decreased  by 
a  factor  of  two  compared  to  the  land  surface  stress  (30%  decrease  of  friction 
velocity,  Figure  3a, b),  and  the  offshore  flow  typically  accelerates  1-2  ms-1 
over  the  fetch  of  a  few  kilometers.  The  exception  is  for  small  stress  over  land 
and  instability  over  the  sea  (Figure  3b),  when  the  stress  is  sometimes  greater 
over  the  sea  than  over  the  land.  Excluding  these  cases,  the  stress  over  the  sea 
is  highly  correlated  with  the  stress  at  the  landmast,  in  spite  of  some  variation 
of  fetch  and  upstream  land  conditions.  However,  much  of  this  correlation  is 
through  the  dependence  of  the  stress  on  wind  speed. 

The  characteristics  of  the  offshore  flow  are  more  systematically  related  to 
travel  time  than  fetch  in  terms  of  the  scatter  in  Figures  7,  8,  12  and  13  in  the 
following  sections,  suggesting  that  the  flow  is  influenced  by  an  internal  decay 
time  scale.  The  turbulence  may  decay  more  near  the  surface  where  the  travel 
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time  to  the  tower  is  longer  (weaker  wind  near  the  surface)  and  the  dissipation 
time  scale  is  shorter  (smaller  turbulent  length  scale  near  the  surface).  For 
offshore  flow,  the  travel  times  at  the  offshore  mast  generally  range  between 
a  few  hundred  seconds  and  about  600  seconds.  The  decay  time  of  convective 
turbulence  can  be  estimated  in  terms  of  the  ratio  of  the  vertical  length  scale 
of  the  turbulence  divided  by  the  velocity  scale  of  the  turbulence.  Nieuwstadt 
and  Brost  (1986)  and  Sorbjan  (1997)  provide  specific  formulations  for  the 
case  where  the  velocity  scale  is  the  free  convection  velocity  and  the  length 
scale  is  the  depth  of  the  convective  mixed  layer.  Applying  such  a  relationship 
to  the  few  RASEX  cases  with  convective  conditions  over  land  predicts  a 
decay  time  scale  on  the  “order”  of  ten  minutes,  which  is  consistent  with  the 
observed  variation  of  turbulence  quantities  with  travel  time  in  Sections  4  and 
5.  In  contrast,  existing  formulations  of  the  surface  stress  in  offshore  flow  are 
formulated  in  terms  of  a  nondimensional  fetch  (Perrie  and  Toulany,  1990; 
Vickers  and  Mahrt,  1997).  More  data  is  required  before  a  new  formulation 
can  be  constructed  in  terms  of  travel  time. 

4.  Upward  buoyancy  flux  from  the  sea  surface 

We  now  analyze  the  flow  of  cool  air  over  a  warmer  surface  in  terms  of  three 
periods  when  the  buoyancy  flux  at  the  surface  exceeds  0.01  °Cms~ 1  for  more 
than  a  half  day  (Table  1).  All  three  cases  axe  characterized  by  very  weak 
instability.  The  time-height  cross-section  (Figure  4)  shows  diurnal  variation 
with  maximum  upward  buoyancy  flux  in  the  morning  due  to  advection  of 
cool  air  from  land.  The  vertical  profiles  are  averaged  over  all  of  the  one-hour 
records  for  each  of  the  three  cases  (Figure  5).  The  turbulence  energy  budget  is 
averaged  over  the  same  data.  The  residuals  of  the  turbulence  energy  budgets 
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(Figure  6)  are  reasonably  small  considering  the  possible  substantial  errors 
in  certain  terms  of  the  budget  and  omission  of  the  pressure  transport  term 
(Section  2.1).  Error  bars  for  the  averaged  profiles  axe  not  shown  because 
much  of  the  variation  within  the  averaging  period  is  due  to  either  diurnal 
trend  or  other  nonstationarity,  rather  than  random  variations.  The  analysis 
in  this  section  also  considers  statistics  based  on  all  of  the  individual  one-hour 
records  within  the  entire  field  program  when  z/L  <  —0.1. 

Table  1.  CASE  STUDIES.  B.F.  is  the  surface  buoyancy  flux  (0(7ms-1), 
“travel”  refers  to  the  travel  time  in  minutes  and  fetch  is  in  kilometers.  Case 
I  for  the  stable  and  unstable  flows  include  time-height  cross-sections  of  the 
buoyancy  flux  (Figures  4  and  9). 


case  study 

Time  (DOY) 

travel 

fetch 

B.F. 

z/L 

unstable-UI 

291.9  -  293.0 

5 

3 

0.02 

10 

unstable-UII 

289.25  -  290 

35 

18 

0.03 

9 

unstable-UIII 

307.7  -  308.35 

4.5 

2.8 

0.01 

10 

stable-SI 

297.0  -  297.7 

5 

2 

-0.02 

7 

stable-SII 

299.0  -  301.0 

6 

2.8 

-0.01 

8 

g|S 

stable-SIII 

303.15  -  303.25 

5 

2 

-0.01 

6.5 

isa 

4.1  Thick  capping  layer  of  downward  buoyancy  flux 

In  RASEX,  offshore  flow  of  cool  air  over  warmer  water  for  short  travel 
times  (near-coastal  zone)  is  generally  characterized  by  a  layer  of  upward 
buoyancy  flux  capped  by  a  relatively  thick  layer  of  downward  buoyancy  flux 
(Case  UI,  Figure  5a  and  Case  UIII,  Figure  5c),  as  also  observed  in  Sun  et 
al.,  (1998)  in  flow  of  cool  air  over  a  warmer  lake.  The  momentum  flux  in 
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the  overlying  layer  of  downward  buoyancy  flux  is  comparable  to  that  near 
the  surface.  These  characteristics  contrast  with  the  textbook  picture  of  a 
convective  internal  boundary  layer  where  the  layer  of  upward  buoyancy  flux 
is  capped  by  a  thin  entrainment  zone  of  downward  buoyancy  flux  and  the 
momentum  flux  decreases  with  height  across  the  convective  boundary  layer. 

The  buoyancy-generation  of  turbulence  energy  is  quite  small  compared 
to  the  shear-generation  and  dissipation  (Figure  6).  The  horizontal  advection 
of  turbulence  energy  is  also  small  in  the  turbulence  energy  budget,  at  least 
near  the  surface  where  it  could  be  evaluated  (Figure  6a,  c).  The  horizon¬ 
tal  advection  term  is  a  little  larger  for  stable  cases  (Figure  11,  Section  5). 
Horizontal  advection  presumably  becomes  increasingly  important  at  higher 
levels  and  also  shoreward  of  the  seamast.  The  main  effect  of  advection  may 
be  through  advection  of  temperature  and  stress.  The  advected  stress  in  turn 
enhances  momentum  transfer  to  the  sea  surface  and  subsequently  enhances 
wave  growth  compared  to  what  would  occur  without  advection  of  stress.  The 
importance  of  advection  of  stress  is  partly  inferred  from  numerical  studies, 
which  will  be  reported  in  a  future  publication. 

For  the  two  short  fetch  cases  (near-coastal  zone  flow),  the  turbulent  trans¬ 
port  term  is  positive  near  the  surface  (Figure  6),  associated  with  convergence 
of  downward  flux  of  turbulence  energy  (not  shown).  For  the  long  fetch  case 
(outer  coastal  zone),  the  vertical  flux  term  is  negative  near  the  surface  corre¬ 
sponding  to  traditional  export  of  turbulence  energy  upward  out  of  the  surface 
layer. 

In  contrast  to  the  short  fetch  cases,  the  flow  in  Case  UII  reaches  the 
tower  after  a  relatively  long  fetch  of  18  km  and  travel  time  of  about  35 
min.  The  offshore  tower  is  now  in  the  outer  coastal  zone  with  respect  to 
the  flow  direction.  The  flow  is  unstable  at  all  levels  and  the  momentum 
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flux  for  Case  UII  (Figure  5b)  is  approximately  constant  with  height,  within 
the  uncertainty  due  to  differences  between  different  sonic  anemometers.  The 
unstable  internal  boundary  layer  is  deep  compared  to  the  tower  layer  because 
of  the  longer  travel  time  over  the  water.  The  buoyancy  flux  is  small  and 
erratic  for  this  case  due  to  small  air-sea  temperature  difference  and  the  flux 
shows  no  obvious  trend  with  height. 

To  further  examine  the  bulk  vertical  structure,  consider  the  ratio  of  the 
buoyancy  flux  in  the  upper  part  of  the  tower  layer  to  that  in  the  lower  part  of 
the  tower  layer  (Eq.  3)  for  all  of  the  records  with  upward  surface  buoyancy 
flux  when  z/L  <  —0.1.  The  buoyancy  flux  ratio  (Figure  7)  for  short  offshore 
travel  times  of  less  than  ten  minutes  (near  coastal  zone)  is  generally  less  than 
one,  indicating  that  the  internal  boundary  layer  is  shallow  and  the  buoyancy 
flux  decreases  significantly  across  the  tower  layer.  In  fact,  the  buoyancy  flux 
ratio  is  negative  for  a  subset  of  the  short  fetch  cases,  as  occurred  in  two  of  the 
case  studies.  With  negative  values,  the  convective  internal  boundary  layer  is 
so  thin  that  the  upper  part  of  the  tower  layer  is  characterized  by  downward 
buoyancy  flux.  The  downward  buoyancy  flux  occurs  over  a  relatively  thick 
layer,  which  probably  extends  above  the  tower  layer. 

For  travel  times  longer  than  10  minutes  (outer  coastal  zone),  the  buoyancy 
ratio  is  always  positive  and  is  typically  of  order  one  with  large  scatter.  That 
is,  the  buoyancy  flux  decreases  only  slowly  with  height  implying  that  the 
internal  boundary  layer  is  deep  compared  to  the  tower  layer. 

The  turbulent  transport  of  turbulence  energy  (Figure  8)  is  often  large  up¬ 
ward  for  unstable  cases  with  short  travel  time  of  less  than  300s  (near  coastal 
zone)  where  advection  from  land  is  most  important.  For  stable  offshore  flow, 
the  pattern  is  quite  different  (next  section). 
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5.  Downward  buoyancy  flux 

Flow  of  warm  air  over  a  cooler  surface  is  now  analyzed  in  terms  of  three 
case  study  periods  where  the  magnitude  of  the  downward  buoyancy  flux 
at  the  surface  is  greater  than  0.01  °Cms~1  for  most  of  the  episode.  The 
time-height  cross-section  for  the  buoyancy  flux  is  shown  for  the  case  with 
the  largest  sustained  downward  buoyancy  flux  (Figure  9),  while  the  vertical 
structure  is  averaged  for  all  of  the  one-hour  records  in  each  of  the  three  cases 
(Figures  10a, b,d,  11  a-c).  All  three  cases  correspond  to  short  fetch  and  short 
travel  time  (Table  1).  We  will  also  analyze  statistics  based  on  all  of  the 
individual  one-hour  records  during  the  entire  field  program  when  z/L  >  0.1. 

For  Case  SI,  the  downward  buoyancy  flux  at  the  sea  surface  is  due  to 
advection  of  warm  air  from  the  heated  land  surface  over  the  cooler  water 
and  exhibits  significant  diurnal  variation.  In  Cases  SI  and  SII,  the  buoyancy 
flux  is  weak  downward  in  the  lowest  20-30  m  and  then  becomes  generally 
negligible  at  higher  levels  (Figures  9,  10a,d). 

5.1  Elevated  stress  maximum 

Above  the  thin  surface  layer  of  downward  buoyancy  flux,  the  stress  and 
vertical  velocity  variance  are  20-60%  larger  than  near  the  surface,  as  can  be 
determined  from  Figure  10  by  squaring  u*  and  <rw.  The  stress  reaches  a  max¬ 
imum  in  the  upper  part  of  the  tower  layer  for  Cases  SI  and  SII  (Figure  10a, c) 
while  the  stress  and  turbulence  energy  slowly  increase  with  height  throughout 
the  tower  layer  for  Case  Sill,  implying  maximum  values  above  the  tower  (Fig¬ 
ure  lOd).  For  shorter  periods,  the  stress  near  the  surface  may  nearly  collapse, 
as  occurred  for  the  one-hour  subperiod  of  Case  SI  (Figure  10b)  where  the 
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measured  stress  and  downward  heat  flux  are  essentially  zero  within  measure¬ 
ment  error.  The  averaged  vertical  structure  in  Cases  SI-SIII  is  quite  different 
from  that  in  traditional  stable  boundary  layers  where  the  stress,  buoyancy 
flux  and  turbulence  energy  decrease  monotonically  with  height  (Caughey  et 
aL,  1979;  Lenschow  et  al.,  1987;  Sorbjan,  1988).  In  the  present  observations 
of  offshore  flow,  the  elevated  maxima  of  stress  and  turbulence  energy  are 
maintained  by  shear-generation  and  presumably  augmented  by  advection  of 
stronger  turbulence  from  land  (Section  5.2). 

The  stress  convergence  below  the  elevated  stress  maximum  acts  to  accel¬ 
erate  the  flow  and  may  account  for  much  of  the  observed  flow  acceleration 
downstream  from  the  coast.  The  observed  vertical  convergence  of  stress  at 
the  sea  mast,  applied  over  the  travel  time  from  the  coast,  corresponds  to  an 
acceleration  of  0.5-1. 0  ms-1.  The  observed  acceleration  ranged  from  0.8  - 
1.3  ms-1.  A  rigorous  assessment  of  the  momentum  budget  is  prevented  by 
inadequate  assessment  of  the  local  horizontal  pressure  gradient. 

The  observed  increase  of  stress  with  height  does  not  appear  to  be  related 
to  instrumentation  differences.  For  long-fetch,  near-neutral  conditions,  the 
observed  stress  is  essentially  constant  with  height  across  the  tower  layer  (or 
decreases  very  slowly  with  height),  as  expected  in  non-ad vective  traditional 
boundary  layers  where  the  boundary  layer  is  much  deeper  than  the  tower 
layer. 

An  elevated  stress  maximum  is  also  observed  by  Glendening  (personal 
communication)  in  an  LES  model  of  flow  from  a  rough  surface  to  a  smooth 
surface  with  zero  buoyancy  flux.  He  used  this  elevated  stress  maximum  to 
define  the  top  of  the  new  internal  boundary  layer.  The  elevated  maximum 
was  maintained  by  horizontal  advection.  Definition  of  the  top  of  the  internal 
boundary  layer  in  terms  of  a  stress  maximum  contrasts  sharply  with  defi- 
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nition  of  the  top  of  the  internal  boundary  layer  in  terms  of  a  minimum  of 
turbulence  cited  in  the  Introduction.  The  differences  of  the  vertical  structure 
may  be  due  to  different  relative  position  downstream  from  the  surface  change. 
The  present  observations  in  offshore  flow  are  short  fetch  (near-coastal  zone) 
and  Mengelkamp  (1991)  finds  that  the  minimum  develops  after  a  certain 
distance  offshore. 

We  now  examine  the  ratio  of  stress  in  the  upper  part  of  the  tower  layer  to 
that  in  the  lower  part  of  the  tower  layer  (Eq.  2)  for  all  of  the  one-hour  records 
where  z/L  >  0.1.  The  momentum  flux  ratio  is  often  large  for  small  travel 
time,  exceeding  two  in  a  significant  fraction  of  the  cases  (Figure  12).  These 
values  correspond  to  a  rapid  increase  of  stress  with  height.  The  momentum 
flux  ratio  decreases  rapidly  with  increasing  travel  time  to  values  closer  to 
unity  for  travel  times  of  ten  minutes  and  longer,  although  the  scatter  is  large. 
For  values  near  unity,  the  stress  decreases  slowly  with  height,  implying  that 
the  internal  boundary  layer  is  deep  compared  to  the  tower  layer. 

5.2  Downward  transport  of  turbulence;  upside-down  boundary 
layer 

For  Cases  SI  and  Sill,  the  vertical  turbulent  transport  of  turbulence  ki¬ 
netic  energy  is  significant  downward  (not  shown),  implying  that  the  main 
source  of  turbulent  kinetic  energy  is  elevated.  This  downward  transport  of 
turbulence  energy  is  consistent  with  the  increase  of  turbulence  and  stress 
with  height,  in  contrast  to  the  usual  boundary  layer.  The  downward  trans¬ 
port  of  turbulence  kinetic  energy  leads  to  significant  vertical  flux  convergence 
of  turbulence  energy  at  the  lower  levels  (Figure  11).  The  turbulence  near  the 
surface  is  therefore  partly  maintained  by  downward  transport  of  turbulence 
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energy.  For  Case  SII,  the  height-dependence  of  the  turbulent  transport  is 
irregular. 

Cases  of  large  elevated  generation  of  turbulence  and  downward  transport 
will  be  referred  to  as  “upside-down”  boundary  layers  in  recognition  of  the 
top-down  transport  of  turbulence  energy.  In  contrast  to  the  classification  in 
Mahrt  (1999),  the  upside-down  boundary  layer  in  the  present  offshore  flow 
occurs  with  relatively  weak  stability  and  is  apparently  partly  maintained 
by  horizontal  advection.  The  downward  transport  of  turbulence  toward  the 
surface  may  be  augmented  by  the  pressure  transport  term  in  the  turbulence 
kinetic  energy  equation  (Smedman  et  al.,  1995). 

Considering  all  of  the  stable  one-hour  records,  the  vertical  transport  of 
turbulence  energy  for  short  travel  times  is  often  large  positive  (Figure  13), 
as  also  Occurred  in  the  unstable  offshore  flow  cases.  These  larger  values 
are  apparently  associated  with  advection  from  land.  In  contrast  to  unsta¬ 
ble  conditions,  the  vertical  transport  of  turbulence  energy  is  sometimes  large 
downward  for  short  travel  times  less  than  300  s  (Figure  13),  again  implying 
that  the  strongest  generation  of  turbulence  is  elevated  in  some  stable  offshore 
flows  in  the  near  coastal  zone.  These  cases  of  downward  transport  of  turbu¬ 
lence  energy  normally  correspond  to  an  increase  of  stress  with  height  (Figure 
14,  momentum  flux  ratio  >  1). 

What  generates  the  upside-down  structure?  Firstly,  the  shear-generation 
term  tends  to  increase  with  height  (Figure  11)  in  contrast  to  the  usual  bound¬ 
ary  layer  where  it  decreases  rapidly  with  height.  Secondly,  the  turbulence 
advected  horizontally  from  land  in  offshore  flow  is  thought  to  decay  more 
slowly  at  higher  levels  where  the  turbulence  length  scale  is  larger  and  travel 
time  is  shorter  (stronger  wind  speed).  The  dissipation  rate  divided  by  the 
turbulence  energy  decreases  with  height  for  both  stable  and  unstable  cases. 
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Finally,  the  advection  of  stress  is  expected  to  increase  with  height.  These 
expectations  are  based  on  numerical  simulations,  to  be  reported  in  a  separate 
manuscript. 

6.  Conclusions 

For  the  present  data,  weakly  convective  internal  boundary  layers  in  flow 
of  cooler  air  over  warmer  water  are  generally  capped  by  a  thick  layer  of 
downward  buoyancy  flux  in  the  near  coastal  zone.  Here,  the  vertically  inte¬ 
grated  buoyancy  flux  is  small  or  even  negative  and  the  vertically-integrated 
turbulence  is  driven  by  shear-generation. 

With  advection  of  warm  air  over  cooler  water  in  the  near  coastal  zone,  the 
vertical  structure  may  be  quite  different  than  that  expected  from  traditional 
conceptual  models  of  the  internal  boundary  layer,  even  though  the  flow  is 
only  modestly  stable  {z/L  <  0.5).  These  cases  most  often  occur  in  the 
afternoon  with  flow  of  warm  air  from  the  heated  land  surface.  In  some  cases, 
the  turbulence  energy  and  stress  increase  with  height,  reaching  an  elevated 
maximum,  and  the  transport  of  turbulence  energy  is  downward  toward  the 
surface.  We  have  referred  to  these  cases  as  “upside-down  boundary  layers”. 
The  increase  of  stress  with  height  appears  to  be  maintained  by  horizontal 
advection  from  land  and  shear-generation  associated  with  accelerating  flow 
over  the  water.  This  upside-down  structure  was  observed  with  travel  times 
less  than  10  minutes  (fetches  generally  less  than  5  km),  here  referred  to  as 
the  near-coastal  zone. 
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Figure  legends 

Figure  1.  Hypothesized  two-layer  structure  for  offshore  flow,  a)  Flow  of 
warmer  air  over  a  cooler  surface.  The  downward-pointing  area  indicates  the 
possibility  of  downward  transport  of  turbulence  energy,  b)  Flow  of  cooler  air 
over  a  warmer  surface,  c)  Small  air-sea  temperature  difference.  The  phase 
speed  of  the  wind  driven  waves  ( Cp )  increases  in  the  offshore  direction.  The 
values  of  offshore  travel  time  for  partitioning  the  near-coastal  zone,  outer 
coastal  zone  and  near-equililbrium  region  are  based  on  analyses  in  Section 
4-5.  Transitions  are  gradual  and  travel  time  values  probably  vary  dramati¬ 
cally  between  different  situations. 

Figure  2.  Relationship  of  the  stability  z/L  to  the  external  Richardson 
number  (Eq.  4)  for  all  of  the  short  fetch  one-hour  records. 

Figure  3.  Relationship  of  the  offshore  friction  velocity  (ms-1)  to  the  land- 
mast  value  for  unstable  (a)  and  stable  (b)  conditions. 

Figure  4.  Time-height  cross-section  for  unstable  case  UI  for  the  buoy¬ 
ancy  flux  (°Cms-1)  where  darker  areas  correspond  to  upward  buoyancy  flux. 
Times  are  GMT  which  is  one  hour  behind  local  solar  time. 
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Figure  5a-c.  Vertical  profiles  of  the  friction  velocity  (ms  x),  buoyancy 
flux  (°Cms-1)  and  aw  (m2s-2)  for  the  three  unstable  cases. 

Figure  6a-c.  The  turbulence  energy  budget  (m2s-3)  for  the  unstable 
cases.  A  is  the  advection  term,  e  the  dissipation,  S  the  shear-generation 
term,  F  the  vertical  divergence  of  the  flux  of  the  turbulence  energy,  B  the 
buoyancy-generation  term  and  R,  the  residual. 

Figure  7.  Buoyancy  flux  ratio  (Eq.  3)  as  a  function  of  travel  time  for  all 
of  the  unstable  records  (zjL  <  —0.1)  . 

Figure  8.  The  dependence  of  the  vertical  flux  of  turbulence  energy  (m3s-3) 
at  18  m  on  travel  time  for  all  of  the  one-hour  records  with  fetch  values  less 
than  5  km  for  unstable  conditions  (z/L  <  —0.1). 

Figure  9.  Time-height  cross-section  of  the  buoyancy  flux  (0Cms-1)  for 
stable  case  SI  where  lighter  areas  correspond  to  stronger  downward  buoyancy 
flux.  Times  are  GMT  which  is  one  hour  behind  local  solar  time. 

Figure  10.  Vertical  profiles  of  the  friction  velocity  (ms-1),  buoyancy  flux 
(°Cms~l)  and  crw  (ms-1)  for  the  three  stable  cases  (a,  c,  d)  and  for  a  one- 
hour  subperiod  during  case  SI  where  fluxes  are  very  weak  (panel  b). 

Figure  11.  The  turbulence  energy  budget  (10-3m2s-3)  for  the  three  stable 
Cases  (a-c).  A  is  the  advection  term,  e  the  dissipation,  S  the  shear-  genera¬ 
tion  term,  F  the  vertical  divergence  of  the  flux  of  the  turbulence  energy,  B 
the  buoyancy-generation  term  and  R,  the  residual. 
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Figure  12.  The  stress  ratio  (Eq.  2)  as  a  function  of  travel  time  for  all  of 
the  one-hour  records  with  stable  conditions  {z/L  >  0.1). 

Figure  13.  The  dependence  of  the  vertical  flux  of  turbulence  energy 
(m3s~3)  at  18  m  on  travel  time  for  all  of  the  one-hour  records  with  sta¬ 
ble  conditions  ( z/L  >  0.1). 

Figure  14.  The  relationship  between  the  turbulence  energy  flux  (m3s-3) 
at  18  m  and  the  momentum  flux  ratio  for  stable  conditions  based  on  bin- 
averaged  values  computed  from  the  one-hour  records.  Also  shown  are  stan¬ 
dard  error  bars. 
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Figure  1.  Hypothesized  two-layer  structure  for  offshore  flow,  a)  Flow  of  warmer 
air  over  a  cooler  surface.  The  downward-pointing  arrow  indicates  the  possibility 
of  downward  transport  of  turbulence  energy,  b)  Flow  of  cooler  air  over  a  warmer 
surface,  c)  Small  air-sea  temperature  difference.  The  phase  speed  of  the  wind- 
driven  waves  (Cp)  increases  in  the  offshore  direction.  The  values  of  offshore 
travel  time  for  partitioning  the  near-coastal  zone,  outer  coastal  zone  and  near¬ 
equilibrium  region  are  based  on  analyses  in  Section  4-5.  Transitions  are  gradual 
and  travel  time  values  probably  vary  dramatically  between  different  situations. 
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Figure  3a.  Relationship  of  the  offshore  friction  velocity  (ms  !)  to  the  land-mast 
value  for  stable  conditions. 
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Figure  4.  Time-height  cross-section  for  unstable  case  UI  for  the  buoy¬ 
ancy  flux  (°Cms-1)  where  darker  areas  correspond  to  upward  buoyancy  flux. 
Times  are  GMT  which  is  one  hour  behind  local  solar  time. 
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Figure  7.  Buoyancy  flux  ratio  (Eq.  3)  as  a  function  of  travel  time  for  all 
of  the  unstable  records  (z/L  <  —0.1)  . 
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Figure  8.  The  dependence  of  the  vertical  flux  of  turbulence  energy  (m3s-3) 
at  18  m  on  travel  time  for  all  of  the  one-hour  records  with  fetch  values  less 
than  5  Wm  for  unstable  conditions  (z/ L  <  —0.1). 
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Figure  9.  Time-height  cross-section  of  the  buoyancy  flux  (°Cms_1)  for 
stable  case  SI  where  lighter  areas  correspond  to  stronger  downward  buoyancy 
flux.  Times  are  GMT  which  is  one  hour  behind  local  solar  time. 
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Figure  10.  Vertical  profiles  of  the  friction  velocity  (ms-1),  buoyancy  flux 
(°Cms~l)  and  aw  (ms-1)  for  the  three  stable  cases  (a,  c,  d)  and  for  a  one- 
hour  subperiod  during  case  SI  where  fluxes  are  very  weak  (panel  b). 
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Figure  12.  The  stress  ratio  (Eq.  2)  as  a  function  of  travel  time  for  all  of 
the  one-hour  records  with  stable  conditions  (z/T  >  0.1). 
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Figure  13.  The  dependence  of  the  vertical  flux  of  turbulence  energy 
(m3s-3)  at  18  m  on  travel  time  for  all  of  the  one-hour  records  with  sta¬ 
ble  conditions  (z/ L>  0.1). 
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Figure  14.  The  relationship  between  the  turbulence  energy  flux  (m3s-3) 
at  18  m  and  the  momentum  flux  ratio  for  stable  conditions  based  on  bin- 
averaged  values  computed  from  the  one-hour  records.  Also  shown  are  stan¬ 
dard  error  bars. 
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10.1  Introduction 

The  coastal  zone  is  of  considerable  practical  importance  since  it  is  a  region  of  intensive 
human  activity  and  rich  in  biological  activity,  compared  to  the  open  ocean.  At  the  same 
time,  the  spatial  variation  of  the  atmospheric  surface  flow  and  the  surface  wave  field  are 
complex  and  both  are  often  in  nonequ ilbrium .  Formulations  of  surface  stress  must 
include  the  influence  of  developing  wave  state  and  shoaling  (Section  3).  Growing  waves 
occur,  with  offshore  flow  and  with  time- varying  wind  fields,  often  induced  by  diurnal 
variation  of  the  heat  flux  discontinuity  at  the  coastal  interface  (Section  5).  Current 
formulations  of  the  surface  flux  may  perform  poorly  in  developing  internal  boundary 
layers  in  offshore  flow  (Section  4).  Difficulties  with  formulation  of  surface  fluxes  in  the 
coastal  zone  is  a  major  topic  of  this  chapter  (Sections  2-3).  Offshore  flow  must  adjust  to 
smaller  surface  roughness  and  different  surface  heat  flux  over  the  sea.  The  modelling 
problem  is  particularly  difficult  with  offshore  flow  of  warm  air  over  a  cool  water 
surface  (Section  4).  The  development  of  stratus  and  other  cloud  regimes  in  the  coastal 
zone  is  beyond  the  scope  of  this  survey.  This  chapter  will  assume  familiarity  with 
previous  reviews  found  in  Donelan  (1990)  and  Geernaert  (1990). 

10.2  Formulation  of  surface  fluxes 

The  surface  stress,  heat  flux  and  moisture  flux  are  formulated  through  the  bulk 
formulas: 


,  >  >•>  —r—t 

(w  u  +w  v 


1/2  ^  _2 
=  Cdu 


( 10-1 ) 


CX.  Geernaert  (ed.).  Air-Sea  Exchange:  Physics,  Chemistry  and  Dynamics,  247-267. 
©  1999  Kluwer  Academic  Publishers.  Printed  in  the  Netherlands. 
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w  0  =  Chu[  0o  —  0] 


(10-2) 


w  q  =  Cg«[ #£>—#] 


(10-3) 


where  Cd  is  the  drag  coefficient,  Ch  is  the  transfer  coefficient  for  heat,  Cq  is  the  transfer 
coefficient  for  moisture,  z  is  the  observational  height,  0,  and  qe  are  surface  aerodynamic 
values  of  the  potential  temperature  and  specific  humidity,  to  be  defined  below,  and  u  is 
the  speed  of  the  vector  averaged  wind  where  the  x-coordinate  has  been  rotated  into  the 
direction  of  the  mean  wind.  Similar  bulk  formulations  can  be  written  for  turbulent 
transfer  of  other  scalars  although  definition  of  the  aerodynamic  surface  quantities 
becomes  problematic  (Sun  et  al.  1998a).  Numerous  variations  of  the  bulk  formula  are 
used  where  the  vector  averaged  speed  is  replaced  by  the  averaged  instantaneous  speed, 
the  total  stress  is  replaced  by  the  stress  in  the  along-wind  direction,  or  the  current 
velocity  is  removed  from  the  wind  (see  Mahrt  et  al.  1998).  Since  the  stress,  wind  and 
wave  propagation  directions  may  be  systematically  different  (Geemaert  1990;  Rieder  et 
al.  1994;  Rieder  1996),  omission  of  the  cross  wind  stress  seems  inadvisable  even  though 
♦his  component  can  be  characterized  by  large  flux  sampling  errors. 

In  the  case  of  weak  winds,  the  precise  definition  of  the  velocity  used  in  the  bulk 
formula  becomes  important  This  velocity  scale  has  been  generalized  to  include  the 
influence  of  "large  convective  eddies"  (Beljaars  1995;  Fairall  et  al.  1996;  Grachev  et  al. 
1998)  or  mesoscale  motions  which  are  on  smaller  scales  than  the  spatial  or 
grid-averaging  scale  (Mahrt  Sun  1995;  Vickers,  Esbensen  1998).  Resolving  this 
problem  is  made  difficult  by  severe  flux  sampling  errors  in  the  weak  wind  case  (Mahrt 
et  al.  1996)  which  can  be  reduced  only  with  long  stationary  records.  With  weak  winds 
and  swell,  the  momentum  flux  may  even  be  transported  from  the  waves  to  the 
atmosphere  (e.g.  Smedman  et  al.  1994). 

The  simplest  formulation  of  the  surface  stress  specifies  the  neutral  drag  coefficient 
to  be  constant  or  a  function  of  wind  speed  but  independent  of  wave  state  (see  studies 
surveyed  in  Geemaert,  1990).  The  transfer  coefficient  for  heat  is  often  specified  to  be 
constant  where  0#  is  equated  to  the  surface  radiation  surface  temperature  of  the  water  or 
the  near-surface  water  temperature  and  qg  is  specified  to  be  the  saturation  specific 
humidity  evaluated  at  the  value  of  the  surface  temperature.  The  appealing  simplicity  of 
this  approach  meets  the  demand  of  economy  required  in  some  large  scale  models. 

However,  in  the  coastal  zone,  stability  strongly  influences  the  transfer  coefficients. 
One  might  expect  stability  -zJL  to  be  smaller  over  water  than  ova:  land  since  surface 
heat  fluxes  are  typically  smaller  over  water  than  over  land,  where  L  is  the  Obukhov 
length 


*z(g/Qv)w  0v 


where  u.  is  the  surface  friction  velocity  (square  root  of  the  surface  stress  magnitude)  and 
k  is  the  von  Kannan  constant  However,  the  shear  stress  is  also  weaker  over  the  sea  than 
over  land  due  to  small  surface  roughness  so  that  the  stability  -zJL  can  significantly 
influence  the  transfer  coefficients  even  without  strong  surface  buoyancy  flux.  As  a 
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result,  stability  may  strongly  influence  the  flux-gradient  relationship  in  the  coastal  zone. 

To  include  the  influence  of  stability  on  the  drag  coefficient  and  transfer  coefficients, 

Monin-Obukhov  similarity  theory  (Monin,  Obukhov  1954)  is  applied  by  either: 

1)  using  similarity  theory  to  reduce  the  drag  and  transfer  coefficients  to  their  equivalent 
neutral  values  (Section  10.2.2)  and  relating  the  neutral  values  to  the  wind  speed  and 
wave  state  or 

2)  explicitly  applying  Monin-Obukhov  similarity  theory  and  relating  the  roughness 
length  to  wave  state  (Section  10.3.2). 


10.2.1  MONIN-OBUKHOV  SIMILARITY  THEORY 

Monin-Obukhov  similarity  theory  provides  a  model  for  the  flux-gradient  relationship  in 
the  surface  layer,  above  the  wave  boundary  layer  (Figure  10.1).  In  the  surface  layer,  the 
flux-gradient  relationship  is  assumed  to  be  independent  of  the  wave  state  and  a 
universal  function  of  only  z/L.  This  flux-gradient  relationship  is  posed  in  terms  of  the 
nondimensional  vertical  gradient 


(ae/aZ)(Kzu*)  (io-5) 

(z/L)  e - — - 

w0(z) 

<| >« ( Z/L)  =  lH>  (  10-6  ) 

u* 

which  have  been  empirically  fitted  to  dependencies  on  z/L  (Hogstrom  1988).  The 
nondimensional  gradient  <j)  can  be  interpreted  as  the  inverse  of  the  mixing  efficiency. 
The  transfer  of  moisture  and  other  scalars  are  often  assumed  to  follow  the  same 
dependence  as  that  for  heat  This  interpretive  survey  will  concentrate  on  transfer  of  heat 
and  momentum  as  examples. 

The  surface  drag  coefficient  and  transfer  coefficient  for  heat  are  estimated  by 

vertically  integrating  the  <{>  functions  downward  to  the  surface  roughness  lengths  in 
which  case  one  can  derive: 


ln(z/zo)  ~  Ym 

C  _  [  _ _ ^ _ -I  r  _ k _ 

h  ln(z/Zo)  ~V/n  ln(z/zoh) 


(10-7) 

(10-8) 


where  za  and  are  the  roughness  lengths  for  momentum  and  heat,  respectively,  and  \j/ 
and  V*  are  the  integrated  analytical  forms  of  the  nondimensional  gradients  for 
momentuni  and  heat,  <j>m  (z/L)  and  «j>4  (z/L).  This  vertical  integration  (Paulson,  1970) 
requires  that  the  fluxes  and  the  wind  direction  are  approximately  independent  of  height 
and  requires  additional  mathematical  approximations  (Enriquez,  Friehe  1997).  The 
nondimensional  vertical  gradients,  <j>M  (z/L)  and  (z/L),  have  not  been  extensively 
evaluated  over  the  sea.  Recently,  Hpjstrup  (1998)  found  that  for  onshore  flow  in  the 
coastal  zone,  the  land-based  stability  dependence  for  4>„  (z/L)  approximated  the  observed 
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values  of  <j>„  although  the  scatter  was  large.  Conflicting  evidence  is  noted  in  the  next 
subsection. 

Note  that  Monin-Obukhov  similarity  theory  does  not  describe  the  actual 
flux-gradient  relationship  below  the  surface  layer  in  the  wave  boundary  layer.  The 
extrapolated  wind  profile  described  by  similarity  theory  vanishes  at  the  roughness 
height  z0  instead  of  at  the  surface.  As  a  result,  the  extrapolated  and  actual  wind  profiles 
in  die  wave  boundary  layer  are  not  expected  to  be  the  same.  The  extrapolated  profile  of 
potential  temperature  is  also  expected  to  be  different  from  the  actual  profile  in  the  wave 
boundary  layer,  and  reaches  the  "aerodynamic"  surface  temperature  at  the  roughness 
height  for  heat,  za-  Both  the  aerodynamic  surface  temperature  0O  and  are  unknown. 
Specifying  0,  to  be  the  surface  radiation  temperature  or  near-surface  water  temperature 
redefines  the  roughness  height  for  heat  za.  Over  the  open  ocean,  the  thermal  roughness 
height  is  normally  related  to  the  roughness  length  for  momentum  (Lui  et  al.  1979).  In 
the  coastal  zone,  Mahrt  et  al.  (1998)  find  little  relationship  between  z*  and  za  where 
zM  appears  to  be  strongly  influenced  by  development  of  internal  boundary  layers  while 
z,  is  more  influenced  by  wave  state.  Both  the  physical  meaning  of  the  roughness  height 
for  heat  and  its  dependence  on  other  parameters  remain  elusive,  especially  in  the  coastal 
zone. 
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Figure  10.1  Idealized  layering  of  the  lower  boundary  layer.  Using  10  ra  as  a  reference  height,  the 
first  scenario  shows  die  ideal  case  where  Monin-Obukhov  similarity  theory  applies.  In 
the  second  example,  the  reference  level  is  in  the  wave  boundary  layer  and  wave  length 
scales  are  required  to  describe  die  local  flux-gradient  relationship.  In  the  third  example, 
the  reference  level  is  above  the  surface  layer  and  bulk  boundary-layer  scaling  is 
required.  In  the  fourth  example,  the  influence  of  the  boundary-layer  depth  extends 
downward  to  the  wave  boundary  layer  and  the  conditions  for  Monin-Obukhov  similarity 
theory  are  not  met  at  any  level. 


10.2.2  REDUCING  TO  NEUTRAL 

To  eliminate  the  influence  of  stability,  the  drag  coefficient  and  transfer  coefficients  are 
sometimes  reduced  to  their  neutral  values  (Deardorff  1968).  This  procedure  attempts  to 
eliminate  the  influence  of  stability  so  that  the  neutral  drag  and  transfer  coefficients  can 
be  studied  as  a  function  of  wave  state  and  wind  speed.  The  reduction  of  the  transfer 
coefficient  to  neutral  conditions  using  Monin-Obukhov  similarity  theory  must  impose 
restrictions  on  the  roughness  length  such  as  Chamock's  relationship  with  constant 
coefficient  (Geemaert,  Katsaros  1986;  Geemaert  1990).  Smith  (1980,  Figure  13), 
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Geemaert  (1988)  and  Mahrt  et  al.  (1996)  find  that  the  reduction  of  the  drag  coefficient 
to  neutral  values  does  not  completely  remove  the  influence  of  stability  although  it  is  not 
clear  if  the  stability  functions  §m{zJL)  and  (j )k(z/L)  themselves  are  incorrect,  measurement 
errors  are  large  or  influences  not  included  in  the  Obukhov  length  are  large. 

10.2.3  WAVE  BOUNDARY  LAYER 

The  wave  boundary  layer  is  the  layer  adjacent  to  the  wave  surface,  but  above  the  very 
thin  laminar  sublayer  at  the  surface  (Figure  10.1).  In  the  wave  boundary  layer,  part  of 
the  atmospheric  transporting  motions  are  directly  induced  by  the  waves  and  associated 
perturbation  pressure  field  in  the  air  (Chalikov,  Belevich  1993;  Hare  et  al.  1997).  The 
depth  of  the  wave  boundary  layer  is  thought  to  scale  with  the  surface  wavelength 
(Chalikov,  Belevich  1993).  Multiple  wave  boundary  layers  associated  with  wind  driven 
waves  and  swell  may  coexist.  However,  the  transport  by  eddies  in  phase  with  the  swell 
may  be  small  since  the  swell  are  characterized  by  small  slopes.  In  the  wave  boundary 
layer,  Monin-Obukhov  similarity  theory  does  not  describe  the  local  flux-gradient 
relationship  which  depends  partly  on  the  amplitude  and  wavelength  of  the  dominant 
surface  wave  (Large  et  al.  1995).  Apparently,  the  profile  functions,  <()„  and  <j)4  depend  on 
zJX  as  well  as  zJL  (or  equivalently  zJL  and  7JL )  where  X  is  a  dominant  wave  dimension, 
either  wave  length  or  wave  height 

Consequently,  Monin-Obukhov  similarity  theory  must  be  applied  in  the  surface 
layer  above  the  wave  boundary  layer,  along  with  specified  roughness  lengths,  in  order 
to  predict  the  surface  fluxes.  Analogously,  over  land,  Monin-Obukhov  similarity  must 
be  applied  in  the  surface  layer  above  the  roughness  sublayer  to  predict  surface  fluxes. 
Eqs.  (10-l)-(10-2)  correctly  predict  surfaces  fluxes  using  Monin-Obukhov  similarity 
theory  for  the  drag  and  transfer  coefficients  (Eqs.  (10-7)-(10-8))  only  if  the  roughness 
lengths  can  be  appropriately  specified. 

There  are  two  fundamental  differences  between  the  flux-gradient  relationship  over  the 
sea  and  that  over  land: 

1)  In  the  roughness  sublayer  over  land,  the  time-averaged  flow  may  vary  horizontally  on 
the  scale  of  the  roughness  elements  due  to  semi-stationary  pressure  perturbations 
anchored  to  individual  roughness  elements.  Over  the  sea,  the  roughness  elements 
(surface  waves)  propagate  so  that  such  microscale  spatial  variability  of  the 
time-averaged  flow  does  not  normally  exist 

2)  Over  the  sea,  the  roughness  length  for  momentum  varies  substantially  with  wave 
state.  Over  land,  z0  is  normally  considered  to  be  independent  of  time  for  a  given  wind 
direction. 

However,  there  is  no  obvious  reason  why  these  differences  would  reduce  the 
applicability  of  Monin-Obukhov  similarity  theory  over  the  sea,  provided  that  the 
appropriate  surface  roughness  lengths  can  be  specified  in  a  reasonably  simple  fashion. 
In  fact  Monin-Obukhov  similarity  theory  is  more  likely  less  applicable  over  land  where 
the  required  assumption  of  homogeneity  is  normally  violated  to  some  degree.  Even  over 
apparently  homogeneous  land  surfaces,  microscale  variations  of  soil  moisture  can 
influence  local  eddy  structure.  Nonetheless,  there  are  special  processes  in  the  coastal 
zone  which  violate  assumptions  required  for  Monin-Obukhov  similarity  theory. 
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10.2.4  THEORY  BREAKDOWN  IN  THE  COASTAL  ZONE 

Existing  similarity  theory  may  break  down  in  the  coastal  zone  due  to  the  following 
influences: 

1.  Strong  horizontal  advection.  leads  to  significant  vertical  divergence  of  the  flux.  For 
example  with  steady-state  temperature  advection,  ude/dx  ■=  -  dw'e/dz-  Then  the 
assumption  of  height-independent  fluxes  is  not  a  good  approximation  and  the  vertical 
integration  of  the  nondimensional  gradients,  <(>,  to  obtain  the  stability  corrections  for 
the  transfer  coefficients,  \|/  (Eqs.  (10-7)-(10-8)),  is  not  valid. 

2.  If  the  flux  decreases  to  small  values  at  the  top  of  the  thin  internal  boundary  layers  in 
offshore  flow,  then  the  vertical  divergence  of  the  flux  is  large.  As  a  result,  standard 
observational  levels,  such  as  10  m,  may  be  above  the  thin  surface  layer  where  one  can 
neglect  the  height-dependence  of  the. flux  (Figure  10.1).  The  depth  of  the  surface 
layer  is  bounded  by  some  small  fraction  of  the  boundary-layer  depth,  sometimes 
chosen  as  10%,  such  that  the  height-dependence  of  the  flux  can  be  neglected  in  the 
surface  layer.  For  example,  in  an  internal  boundary  layer  of  50  m  depth,  the  surface 
layer  would  theoretically  be  less  than  5  m  depth.  Note  that  influences  1)  and  2)  are 
not  independent 

3.  With  offshore  flow  of  warm  air  over  cold  water,  the  turbulence  may  be  a  top-down 
process  where  the  primary  source  of  turbulence  is  above  the  surface  inversion 
(Section  10.4.2).  It  is  not  dear  if  Monih-Obukhov  similarity  theory  is  valid  in  this 
"upside-down  boundary  layer". 

4.  The  turbulence  in  offshore  flow  does  not  establish  equilibrium  with  the  rapidly 
evolving  mean  flow  immediately  downstream  from  the  coast. 

5.  With  cold  (warm)  air  advection  in  offshore  flow,  the  wind  vector  tends  to  rotate  to  the 
right  (left)  with  height  and  the  surface  stress  vector  is  directed  to  the  right  (left)  of  the 
surface  wind  vector  (Geemaert  1988).  Integration  of  Monin-Obukhov  similarity 
theory  assumes  that  the  wind  and  stress  vectors  are  aligned  and  their 
height-dependence  can  be  neglected. 

6.  It  may  be.  that  Monin-Obukhov  similarity  theory  is  valid  but  the  stability  functions,  ^ 
and  §k,  are  not  correctly  calibrated. 

If  the  internal  boundary  layer  is  sufficiently  thin,  yet  the  wave  boundary  layer  is  deep 
(large  surface  wavelength),  the  surface  layer  may  be  "squeezed  out"  as  postulated  by 
Mahrt  et  al.  (1998)  and  shown  in  Figure  10.1  (right  hand  side).  That  is,  there  is  no  layer 
where  the  flux-gradient  relationship  depends  only  on  zJL  and  the  influence  of 
boundary-layer  depth,  h,  and  surface  wave  dimension  X  must  be  included  as  additional 
scaling  variables.  Grant  (1992)  suggests  that  ty(z/L)  for  the  near  neutral  boundary  layer 
should  be  generalized  to  be  of  the  form  <|)(z//t,  h/L ,  ujfh)  where  /is  the  Coriolis 
parameter.  Khanna  and  Brasseur  (1997)  consider  the  form  <J»(z/L,  h/L).  In  the  LES 
results  of  Khanna  and  Brasseur  (1997),  the  nondimensional  shear  decreases  above  the 
surface  layer,  as  is  also  observed  by  Smedman  and  Johannson  (1997)  in  shallow 
offshore  boundary  layers.  With  a  well-mixed  interior  of  the  boundary  layer,  the  vertical 
gradient  decreases  faster  with  height  than  the  flux.  Mahrt  et  al.  (1998)  find  that  the  thin 
depth  of  observed  offshore  internal  boundary  layers  suppresses  heat  transfer  by  the  large 
convective  eddies  which  in  turn  implies  that  <|>Ais  larger  than  predicted  by 
Monin-Obukhov  similarity  theory  and  depends  on  z/h. 
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Davidson  (1974)  found  that  and  the  drag  coefficient  are  both  a  function  of 
stability  and  wave  age  and  that  these  dependencies  are  difficult  to  sort  out  since  the 
wave  age  and  stability  were  significantly  correlated  for  their  data.  Bergstrom  and 
Smedman  (1995)  examined  the  functional  dependence  ^  (z/L,  CJu.)  but  found  that  the 
relationship  to  wave  state  was  statistically  insignificant,  although  they  note  that  their 
data  represents  a  relatively  narrow  range  of  conditions.  Here,  C  is  the  phase  speed  of 
the  dominant  wave.  ' 

Since  for  offshore  flow,  there  may  be  no  level  where  <j>w  is  a  function  of  only  z/L,  it  is 
useful  to  consider  a  more  general  formulation  such  as 

*  (3muZ). ( —  Z/X,  Cplu,,zlh)  (10-9) 

Note  that  the  arguments  z/L,  zfk,  CJu.,  zJh  are  not  independent.  Large  et  al.  (1995) 
prefer  to  partition  the  nondimensional  shear  into  two  functions,  the  traditional 
(z/L)  and  a  modifying  function  %(zfk),  as  also  pursued  by  Vickers  and  Mahrt  (1998)" 
Specific  forms  of  Eq.  (10-9)  based  on  actual  eddy  correlation  data  have  not  been 
established. 

Monin-Obukhov  similarity  theory  is  also  complicated  by  the  possible  dependence  of 
the  von  Karmen  "constant"  on  the  roughness  Reynolds  number.  Re.  defined  as  u.  zjv. 
Generally,  k (Re*)  is  found  to  decrease  with  increasing  roughness  Reynolds  number 
(Tennekes,  1968).  With  this  possibility,  evaluation  of  $m(zJL)  from  data  using  Eq.  (10-6) 
results  in  two  unknowns,  k (Re.)  and  §m(zJL).  Oncley  et  al.  (1996),  Miller  et  al.  (1997) 
and  others  eliminate  <t>„(z/L)  as  an  unknown  by  considering  near  neutral  cases  and 
assuming  4>m(z/L)  =1.  Then,  x(Re.)  is  estimated  as 

K (Re*)  (10-10) 

zdu/dz 


Since  both  Re.  and  z/L  depend  on  their  dependencies  are  difficult  to  isolate  from 
data.  For  example,  failure  to  include  the  dependence  k (Re.),  if  important,  would  alter 
the  value  of  §m(z/L)  computed  from  data.  Conversely,  if  influences  of  z/L  are  not 
completely  negligible,  such  influences  could  create  an  artificial  dependence  of  Kon  the 
roughness  Reynolds  number.  The  problem  becomes  even  more  complex  if  the 
nondimensional  shear  is  dependent  on  wave  age  (Eq.  (10-9)  since  the  wave  age  also 
depends  on  u.. 

103  Wave  state 

In  the  coastal  zone,  the  stress  is  influenced  by  shoaling  processes  and  wave  breaking  as 
waves  propagate  into  shallow  water  (Smith  1980;  Freilich  et  al.  1990).  For  waves 
propagating  toward  the  shore,  wave  modification  first  occurs  when  the  depth  of  the 
water  shallows  to  be  about  1/4  of  the  dominate  wavelength.  At  this  point,  the  wave 
amplitude  begins  to  increase,  the  wavelength  and  propagation  speed  begin  to  decrease 
and  the  spectral  characteristics  change.  Closer  to  the  shore  where  the  depth  becomes 
less  than  one  wave  height,  dramatic  wave  steepening  and  breaking  occur  (Thornton, 
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Guza  1982,  1983;  Holman,  Sallenger  1985;  Holland  et  al.  1995).  Irregularities  of  the 
bottom  topography  along  the  coast  and  wave  refraction  lead  to  irregularities  in  the 
surface  wave  field  along  the  coast  (Munk,  Traylor  1947). 

Information  on  wave  state  is  necessary  for  modelling  surface  stress  in  the  coastal 
zone.  This  information  is  included  at  three  levels  of  approximation:  1)  relating  the  drag 
coefficient  or  roughness  length  to  the  wave  age  is  the  most  common  approach  (Section 
10.3.1),  2)  relating  the  stress  to  different  frequency  bands  of  the  wave  field,  such  as 
swell  and  wind-driven  modes,  is  more  complete  but  requires  more  information  on  the 
wave  field  (Section  10.3.2)  and  3)  modelling  the  wave  age  in  terms  of  fetch  is  less 
accurate  but  is  useful  when  explicit  information  on  wave  state  is  not  available  (Section 
10.3.3). 

10.3.1  WAVE  AGE  DEPENDENCE 

With  offshore  flow,  the  wind  driven  waves  are  growing  in  the  downstream  direction 
which  become  evident,  seaward  of  the  inner  shoaling  zone.  A  number  of  investigators 
have  documented  that  the  stress  is  greater  over  a  young  and  developing  wave  field  than 
over  an  older  wave  field,  which  is  in  near-equilibrium  with  the  wind  field  (e.g.  Donelan 
1990;  Kitaigorodskii  1973;  Nordeng  1991;  Geemaert  et  al.  1987,  1988;  Donelan  1990; 
Maat  et  al  1991).  At  least  two  mechanisms  contribute  to  the  dependence  of  stress  on 
wave  age:  1)  younger  waves  propagate  with  slower  phase  speed  relative  to  the  wind  and 
thus  provide  greater  bulk  shear,  and  2)  younger  growing  waves  may  be  steeper,  which 
can  lead  to  enhanced  flow  separation  from  individual  wave  crests.  Younger  developing 
waves  occur  with  changing  wind  vector  as  well  as  with  fetch  limited  off-shore  flow. 

Previous  studies  suggest  that  the  drag  coefficient  is  a  maximum  for  a  wave  age  of  10 
(Nordeng  1991)  and  7  (Geemaert,  Smith  1996).  Kitaigorodskii  et  al.  (1995)  similarly 
find  maximum  roughness  length  at  an  intermediate  wave  age.  In  the  absence  of  swell 
(inland  seas),  the  drag  coefficient  may  be  immeasurably  small  immediately  downstream 
from  the  shore  where  waves  have  not  yet  developed  or  are  very  small  amplitude.  The 
subsequent  wave  growth  in  the  downstream  direction  leads  to  larger  stress,  as  sketched 
in  Figure  10.2.  The  wave  growth  in  the  downwind  direction  can  be  further  enhanced  by 
acceleration  of  the  offshore  wind  resulting  from  the  smaller  surface  roughness  over  the 
sea  than  over  land.  However,  at  some  point,  the  increasing  phase  speed  of  the  waves 
and  resulting  reduction  of  the  relative  flow  of  the  wind  over  the  waves  becomes  more 
important  than  the  effect  of  increasing  wave  amplitude  and  any  increase  of  slope.  Then 
the  stress  and  drag  coefficient  begin  to  decrease  further  downstream.  Consequently,  the 
drag  coefficient  reaches  a  maximum  at  an  intermediate  wave  age  and  decreases  with 
further  increase  of  wave  age.  Similar  effects  appear  in  the  model  of  Hansen  and  Larsen 
(1997)  where  the  Chamock  constant  reaches  a  maximum  at  a  wave  age  of  about  5.  Most 
observations  occur  at  wave  ages  greater  than  this  intermediate  wave  age  so  that  the  drag 
coefficient  is  generally  considered  to  decrease  with  increasing  wave  age,  as  in  the 
studies  surveyed  below. 

Vickers  and  Mahrt  (1997)  find  that  onshore  flow  with  shoaling,  can  also  occur  with 
small  values  of  the  wave  age.  The  shoaling  reduces  the  wave  phase  speed  and  decreases 
the  numerical  value  of  the  wave  age.  Consequently,  shoaling  complicates  the  physical 
interpretation  of  wave  age  in  the  coastal  zone.  The  corresponding  drag  coefficient  may 
increase  by  more  than  50%  with  such  shoaling.  As  a  result,  Kitaigorodskii  et  al.  (1995) 
related  the  momentum  roughness  length  to  wave  breaking  characteristics.  Mahrt  et  al. 
(1998)  found  no  detectable  increase  in  the  transfer  coefficients  for  heat  and  moisture 
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with  wave  breaking.  Relative,  insensitivity  of  the  transfer  coefficients  for  heat  and 
moisture  to  wave  state  and  wind  speed  was  also  found  in  Smith  (1980),  Makin  and 
Mastenbroek  (1996)  and  references  therein. 

Since  the  stress  appears  to  be  related  to  the  relative  flow  of  the  air  with  respect  to  the 
phase  speed  of  the  dominate  waves  (Kitaigorodskii  1973),  one  is  tempted  to  re-define 
the  drag  coefficient  in  terms  of  the  relative  flow  u — Cp  where  Cp  is  ideally  the 

component  of  the  phase  speed  of  the  dominant  wave  in  the  wind  direction.  However, 
the  stress  does  not  necessarily  vanish  as  u  —  Cp  approaches  zero,  since  waves  with 

frequencies  different  from  the  dominant  wave  contribute  to  the  stress  and  propagate 
with  different  phase  velocities.  For  example,  with  mature  waves,  much  of  the  stress  is 
thought  to  be  associated  with  capillary  waves  and  therefore  not  related  to  the  phase 
speed  of  the  dominant  waves.  As  a  result  of  the  stress  from  capillary  waves,  the 
alternate  drag  coefficient,  defined  as  the  ratio  of  the  stress  to  u-cp,  approaches  infinity 

as  u- Cp  approaches  zero. 

However  the  traditional  drag  coefficient  computed  from  (Eq.  (10.1))  may  be  better 
related  to  u-  Cp  than  u  alone.  To  examine  this  problem,  we  computed  the  drag 

coefficient  from  eddy  correlation  data  collected  at  a  tower  2  km  off  the  Danish  coast  in 
the  Ris0  Air  Sea  Experiment  (RASEX,  Barthelmie  1994;  Hpjstrup  et  al.  1994;  Mahrt  et 
al.  1996).  The  drag  coefficient  reaches  a  minimum  value  near  u-Cp=  0  m/s  and  a 

regression  model  of  the  drag  coefficient  based  on  u- Cp explains  more  variance  than 

that  based  on  u.  However,  the  drag  coefficient  is  more  closely  related  to  wave  age  than 

«-  Cp  probably  because  wave  age  implicitly  accounts  for  the  relative  flow  over  the 

waves  as  well  as  implicitly  includes  the  effect  of  wave  amplitude  and  steepness. 

Geemaert  et  al.  (1987)  proposed  a  model  of  the  dependence  of  the  neutral  drag 
coefficient  on  wave  age  of  the  form 

Cdn=  b(Cp/u*)~2'3  (10-11) 

where  b  is  on  the  order  of  10'2.  Vickers  and  Mahrt  (1997)  found  that  this  model 
performed  well  in  the  coastal  zone  even  after  accounting  for  the  role  of  self-correlation. 
They  found  that  -2/3  was  indeed  the  best  fit  to  the  exponent  while  the  best  fit  value  of 
b  was  approximately  7  xlO  s. 

The.  drag  coefficient  can  also  be  computed  directly  from  Monin-Obukhov  similarity 
theory  in  which  case  the  roughness  length  is  a  function  of  wave  state.  The  Chamock 
(1955)  prediction  of  the  roughness  length,  0.019  (  u.2  /  g),  is  often  generalized  to  include 
a  dependence  on  inverse  wave  age  (Toba,  Koga  1986;  Maat  et  al.  1991;  Donelan  1990; 
Smith  et  al.  1992).  If  we  also  incorporate  the  smooth  flow  contribution  to  the  roughness 
length  (Donelan  1990;  Fairall  et  al.  1996),  then  the  full  relationship  is  of  the  form 

Zo=K(ul  / g )  (u*/Cp)P+0.UMu*  (10-12) 
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Figure  10.2  Plausible  schematic  evolution  of  the  internal  boundaiy  layer  for.  unstable  and  very  stable 
offshore  flows.  In  the  unstable  case,  the  x-axis  can  also  be  qualitatively  interpreted  as  the 
downstream  distance. 


where  typically  K  =0.48  and  p  —1  (Smith  et  al.  1992),  and  v  is  the  viscosity  of  air.  Wu 
(1968)  suggested  an  additional  term  due  to  parameterized  surface  tension  which  was 
recently  applied  in  Alam  and  Curry  (1997). 

Alternatively,  the  model  of  Kitaigorodskii  (1970)  is  of  the  form, 

Zo=Kcs&x p(  ~K.Cpf  u #)  (10-13) 

where  K  =  0.3,  <r  is  the  root  mean  square  amplitude  of  the  waves  and  k  is  von  Kannan's 
constant.  The  Kitaigorodskii  representation  follows  from  explicitly  including  the  shear 
between  the  wind  and  the  wave  velocity  and  integrating  over  the  full  wave  spectrum. 

A  significant  fraction  of  the  variance  explained  by  these  models  is  due  to 
self-correlation  since  the  drag  coefficient  and  roughness  length  are  both  defined  in 
terms  of  the  surface  friction  velocity  and  are  both  formulated  in  terms  of  the  surface 
friction  velocity  on  the  right  hand  sides  of  Eqs.  (10-11)-(10-13)  (e.g.  Smith  et  al.  1992; 
Vickers,  Mahrt  1997).  To  avoid  self-correlation,  the  wave  age  is  sometimes  expressed 
in  terms  of  Cpfu  where  u  is  the  wind  speed  at  a  standard  level  or  at  a  fixed  height 

relative  to  the  surface  wavelength  (Donelan  1990). 

The  roughness  length  is  sometimes  related  to  other  characteristics  of  the  wave  field. 
For  example,  Anctil  and  Donelan  (1996)  relate  the  spatial  variation  of  the  roughness 
length  in  the  shoaling  zone  to  the  wave  age,  root  mean  squared  displacement  height  of 
the  waves  and  root  mean  square  wave  slope.  Kitaigorodskii  et  al.  (1995)  relate  the 
roughness  length  to  the  effective  wave  height. 

10.3.2  MULTI-MODE  MODELS 

The  need  for  models  based  on  multiple  wave  modes  is  motivated  by  the  fact  that  the 
wind  driven  waves  and  swell  often  propagate  in  different  directions.  The  swell  modifies 
the  stress  direction  so  that  it  may  be  different  from  the  wind  direction  (Geemaert  et  al. 
1995).  A  family  of  parameterizations  of  the  roughness  length  as  a  weighted  integral  of 
the  wave  spectra  have  been  formed  from  the  original  framework  of  Kitaigorodskii 
(1973)  as  in  Hansen  and  Larsen  (1997).  See  also  Chalikov  and  Belevich  (1993)  and 
papers  surveyed  in  Geemaert  (1990).  Kitaigorodskii  et  al.  (1995)  focus  on  the  width  of 
the  wave  dissipation  regime  in  an  effort  to  understand  the  variability  of  surface 
roughness.  The  most  complete  description  of  the  wave  field  can  be  obtained  from 
explicit  wave  models  such  as  the  "WAM"  model  (WAMDI  group,  1988).  The 
advantage  of  relating  the  surface  roughness  to  the  wave  spectra  avoids  the 
parameterized  coupling  between  the  surface  roughness  and  with  the  wind  field  or  stress 
field  itself  that  occurs  with  the  approaches  in  the  preceding  section.  At  the  same  time, 
complete  information  on  the  full  wave  spectra  is  normally  not  available. 

Simplified  models  which  allow  some  information  of  differences  between  wave 
modes  include  the  two  band  model  of  Donelan  (1982)  which  distinguishes  between  the 
drag  associated  with  short  and  long  waves.  Vickers  and  Mahrt  (1997)  develop  a  crude 
indicator  of  the  wave  state  based  on  that  band  width  of  the  spectra  which  accounts  for 
50%  of  the  wave  energy.  Narrow  band  spectra  are  associated  with  near-equilibrium 
single  mode  waves  and  smaller  drag  coefficient  while  broad  band  spectra  are  more 
associated  with  multi-peaked  spectra,  confused  seas  and/or  non-equilibrium  wave  state, 
all  leading  to  larger  drag  coefficients. 
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Here  X  is  the  upwind  fetch  distance  and  wc  and  a  are  determined  empirically.  Eq.  (10- 
14)  explains  more  variance  than  Eq.  (10-15)  but  is  characterized  by  self-correlation. 
Vickers  and  Mahrt  (1997)  find  that  Eq.  (10-15)  is  a  reasonable  approximation  for 
offshore  flow  in  the  coastal  zone  flow,  although  errors  are  larger  compared  to 
formulations  based  on  wave  state  information. 

Direct  relationship  of  the  roughness  length  to  the  fetch  can  be  found  in  Geemaert 
(1988b).  More  complete  models  first  describe  the  dependence  of  the  wave  field  on  fetch 
and  then  link  the  drag  to  the  fetch-dependent  wave  field  (Geemaert  1990). 

10.4  Internal  boundary  layers 

With  any  wind  component  perpendicular  to  the  coast,  internal  boundary  layers  form  due 
to  the  temperature  and  surface  roughness  contrast  between  the  water  and  land  surfaces. 
The  internal  boundary  layer  is  the  layer  of  air  adjacent  to  the  surface  which  is 
influenced  by  the  new  surface  (e.g.  Garratt  1990). 

10.4.1  UNSTABLE  CASE 

With  offshore  flow  of  cold  air  over  warmer  water,  a  convectively  driven  internal 
boundary  layer  forms  in  response  to  the  upward  buoyancy  flux  over  the  water.  This 
layer  thickens  in  the  downstream  direction  (Figure  10.2).  These  internal  boundary  layers 
are  sometimes  found  to  be  well  defined  by  a  sharp  change  of  properties  between  the 
convectively  driven  turbulence  and  fluid,  at  least  in  terms  of  instantaneous  observations 
(Raynor  et  al.  1979;  Sun  et  al.  1998b).  It  is  not  known  if  such  sharply  defined  tops  of 
the  internal  boundary  layer  are  normal  nor  is  it  known  if  the  internal  boundary  layer  top 
for  time-averaged  flow  is  thick  due  to  vacillation  of  the  boundary-layer  top. 
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Figure  10.3  shows  the  averaged  vertical  structure  for  a  nine-hour  period  of 
stationary  offshore  advection  of  cold  air  observed  in  RASEX  (Section  10.3.1).  Several 
features  occur  in  Figure  10.3  which  are  not  typical  of  the  "textbook”  unstable  internal 
boundary  layer: 

1)  For  the  averaged  flow,  the  entrainment  zone"  with  downward  heat  flux  is  thicker 
than  the  convective  layer  of  upward  buoyancy  flux. 

2)  The  total  buoyancy-destruction  of  turbulence  in  the  entrainment  zone  is  greater  than 
the  buoyancy-generation  of  turbulence  in  the  surface-based  convective  layer. 
Therefore,  turbulence  near  the  top  of  the  internal  boundary  layer  must  be  maintained 
by  local  shear-generation,  as  observed  in  Sun  et  al.  (1998b). 

3)  The  turbulence  and  stress  are  a  maximum  near  the  top  of  the  convective  layer. 

4)  The  mean  shear  increases  with  height  corresponding  to  convex  curvature.  This 
curvature  implies  inflection  points  near  the  surface  and  somewhere  above  the  45  m 
tower  layer. 


Unstable  Internal  Boundary  Layer 


r.  9  <c) 

Figure  10.3  The  5-hour  averaged  vertical  structure  of  relative  stationary  offshore  flow  of  cool  air  over 

warmer  water  observed  in  RASEX.  Fluxes  are  based  on  perturbations  from  a  10-minute 
simple  mean. 
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Is  this  structure  unique  to  these  data  or  does  the  traditional  concept  of  internal  boundary 
layers  not  apply  immediately  downstream  from  the  coast?  The  RASEX  data  set  is  the 
first  data  set  with  detailed  vertical  resolution  of  the  time-averaged  structure  in  offshore 
flow. 

Kallstrand  and  Smedman  (1997)  compare  various  models  of  growth  of  the  internal 
boundary  layer  against  aircraft  data  collected  in  a  developing  internal  boundary  layer 
over  land  with  onshore  flow.  Although  most  of  the  models  have  been  developed  for  the 
unstable  internal  boundary  layer  over  land  (e.g.,  Melas,  Kambezidis  1992;  Gryning, 
Batchevarova  1990),  they  should  in  principal  apply  over  the  sea  in  the  coastal  zone.  The 
internal  boundary  layer  is  traditionally  modelled  on  two  separate  scales  (Garratt  1990): 
small  scales  on  die  order  of  a  kilometer  or  less  and  the  mesoscale  on  a  horizontal  scale 
of  tens  of  kilometers  or  more.  On  the  scale  of  a  kilometer  or  less,  the  heated  internal 
boundary  layer  entrains  upward  through  the  old  boundary  layer  (Vugts,  Businger  1977). 
On  this  scale,  the  initial  growth  rate  of  the  convective  internal  boundary  layer  is  thought 
to  be  linearly  proportional  to  the  strength  of  the  turbulence  in  the  convective  internal 
boundary  layer  such  that  (Hpstrup  1981;  Brutsaert  1982;  Garratt  1990) 

dhIBL  Gw  (10-18 ) 

— =ca 

where  is  the  depth  of  the  internal  boundary  layer  and  a„  is  the  standard  deviation  of 
the  turbulent  fluctuations  of  vertical  velocity,  often  parameterized  in  terms  of  the 
surface  friction  velocity  and  the  convective  velocity  scale.  The  wind  speed  U  is 
evaluated  at  the  top  of  the  internal  boundary  layer.  Since  stratification  of  the  overlying 
fluid  is  neglected,  no  additional  velocity  scales  are  included.  The  value  of  C  is  thought 
to  be  order  of  unity  (Mahrt  1996). 

Some  investigators  have  formulated  analytical  expressions  for  the  small  scale 
growth  of  the  internal  boundary-layer  depth.  Fbr  example,  Andreas  et  al.  (1979,  1981) 
model  the  depth  in  terms  of  the  fetch.  Observations  of  such  internal  boundary-layer 
growth  are  summarized  in  Andreas  and  Cash  (1998). 

Further  downstream,  on  the  mesoscale,  the  convective  internal  boundary  layer 
completes  its  growth  through  the  old  boundary  layer  (Figure  10.2)  and  begins  entraining 
nonturbulent  fluid,  as  occurs  in  flow  of  cold  air  over  warm  water  studied  in  Chang  and 
Braham  (1991),  Smith  and  MacPherson  (1987),  Rogers  et  al.  (1995)  and  Briimmer 
(1996).  A  well-defined  capping  inversion  usually  develops  and  information  on  the 
stratification  of  the  overlying  fluid  is  required  to  predict  further  growth  of  the 
convective  internal  boundary  layer.  Models  for  this  case  are  surveyed  in  Garratt  (1990) 
and  Kallstrand  and  Smedman  (1997). 

10.4.2  STABLE  CASE 

The  growth  of  the  stable  internal  boundary  layer  due  to  flow  of  warm  air  over  cooler 
water  is  expected  to  be  slower  because  of  buoyancy  destruction  of  turbulence  associated 
with  the  downward  heat  flux.  For  example,  based  on  temperature  profiles  in  Gryning 
(1985),  the  growth  rate  dhmJdx  in  the  stable  internal  boundary  layer  in  the  coastal  zone 
is  less  than  1%.  A  thick  residual  layer  of  decaying  turbulence  extends  from  the  top  of 
the  thin  stable  internal  boundary  layer  to  the  top  of  the  advected  continental  boundary 
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layer  (Rogers  et  al.  1995).  The  growth  of  the  internal  boundary  layer  for  stable 
conditions  is  also  described  by  Mulheam  (1981)  and  Garratt  (1990). 

The  vertical  structure  of  the  stable  internal  boundary  layer  may  be  well  defined  in 
offshore  flow  (e.g.  Garratt,  Ryan  1989)  where  turbulence  quantities  monotonically 
decrease  with  height  However,  if  the  water  temperature  is  much  cooler  than  the 
advected  air,  the  turbulence  at  the  surface  may  collapse  resulting  in  decoupling  of  the 
overlying  advected  turbulence  from  the  surface  (Figure  10.2).  Such  collapse  cannot  be 
predicted  with  existing  similarity  theory  which  is  the  basis  for  surface  fluxes  in 
numerical  models.  Above  the  surface  inversion  over  the  water,  the  advected  turbulence 
from  land  decays  in  the  downstream  distance.  Since  the  change  of  surface  heat  flux, 
viewed  by  a  moving  Lagrangian  column,  is  almost  instantaneous  as  it  crosses  the 
land-sea  boundary,  the  forcing  time  scale  due  to  the  change  of  surface  flux  is  small 
compared  to  the  internal  time  scale  of  the  turbulence  (turbulent  length  scale/turbulent 
velocity  scale).  Then  the  turbulence  decays  in  a  self-similar  fashion  leading  to  a  square 
root  dependence  on  time  (Sorbjan  1997). 

In  addition  to  the  decaying  turbulence  above  the  surface  inversion,  elevated 
shear-generation  of  turbulence  is  observed  further  downstream.  The  flow  above  the 
surface  inversion,  which  was  part  of  the  boundary  layer  over  land,  accelerates  and  forms 
a  low  level  jet  (Smedman  et  al.  1995),  analogous  to  formation  of  a  low-level  jet  above 
the  nocturnal  surface  inversion  (Figure  10.2).  Further  downstream,  the  shear  on  the 
underside  of  the  jet  eventually  generates  turbulence  and  re-establishes  a  surface-based 
boundary  layer.  Still  further  downstream,  the  flow  may  become  near  neutral  as  the  air 
finally  cools  to  the  value  of  the  sea  surface  temperature  (Smedman  et  al.  1997). 

The  RASEX  data  for  offshore  flow  of  warm  air  over  cooler  water  (Figure  10.4)  also 
indicate  that  the  classical  concept  of  an  internal  boundary  layer  does  not  apply.  The 
stress,  heat  and  turbulence  energy  are  smaller  near  the  surface  and  reach  maximum 
values  at  higher  levels  in  the  tower  layer.  Note  that  the  strength  of  the  turbulence  for  the 
stable  case  is  comparable  to  that  for  the  unstable  case  (Figure  10.3)  because  the  wind 
speeds  for  the  stable  case  are  larger.  For  shorter  periods,  the  stress  and  turbulence  in  the 
stable  case  may  temporarily  collapse  near  the  surface  in  that  the  stress  is  zero  within 
measurement  error.  Then  the  elevated  turbulence  is  semi-detached  from  the  surface 
corresponding  to  an  upside-down  boundary  layer. 

10.4.3  ROUGHNESS  CHANGE 

Internal  boundary  layers  may  also  form  due  to  spatial  variation  of  the  surface  roughness 
between  the  land,  shoaling  zone  and  open  ocean.  With  offshore  flow,  the  smaller 
roughness  over  the  sea  compared  to  over  the  land  leads  to  flow  acceleration  and 
decreased  turbulence  strength  in  the  downstream  direction,  as  observed  in  Smith  and 
MacPherson  (1987).  Although  the  roughness  increases  in  the  downstream  direction 
(Figure  10.2),  it  remains  much  smaller  than  the  roughness  over  land. 

Normally,  at  least  some  surface  temperature  discontinuity  occurs  at  the  coast  so  that 
internal  boundary  layers  result  as  combination  of  surface  roughness  and  surface  heat 
flux  changes.  For  example,  Barthelmie  et  al.  (1996)  find  that  in  unstable  boundary 
layers  associated  with  flow  of  cool  nocturnal  air  over  warmer  water,  winds  accelerate 
due  to  smaller  roughness  over  the  water.  However  with  stable  internal  boundary  layers 
associated  with  daytime  flow  of  warmer  air  over  the  water,  the  reduction  of  downward 
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Figure  10.4  The  9-hour  averaged  vertical  structure  of  relative  stationary  offshore  flow  of  warm  air  over 
cooler  water  observed  in  RASEX.  Fluxes  are  based  on  perturbations  from  a  10  minute 
simple  mean. 

mixing  of  momentum  due  to  stable  stratification  counteracts  the  effect  of  decreased 
surface  roughness  on  the  wind  speed.  As  a  result,  acceleration  over  the  water  was  not 
normally  observed  in  the  daytime. 

10.5  Local  circulations 

10.5. 1  SEA/LAND  BREEZES 

The  land-sea  interface  not  only  leads  to  development  of  internal  boundary  layers  but 
also  modifies  the  low  level  pressure  field  through  differential  surface  heating.  On  sunny 
days,  surface  heat  flux  from  the  heated  land  surface  warms  the  atmosphere  and  lowers 
the  surface  hydrostatic  pressure  over  the  land.  Simultaneous  heating  does  not  occur  over 
the  water  surface  since:  a)  the  heat  capacity  of  the  water  is  large,  b)  the  solar  radiation  is 
absorbed  over  the  first  few  meters  instead  of  at  the  surface,  as  occurs  over  land,  and  c) 
heat  from  absorbed  solar  radiation  is  mixed  downward  by  oceanic  turbulence  which  is 
much  more  efficient  than  thermal  conduction  in  the  soil.  As  a  result,  a  horizontal 
pressure  gradient  force  develops  directed  toward  the  land  surface.  In  the  absence  of 
significant  opposing  horizontal  pressure  gradient  on  the  synoptic  scale,  the  local 
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horizontal  pressure  gradient  drives  flow  inland  over  the  heated  land  surface,  referred  to 
as  the  sea  breeze  (Pielke  1984;  Simpson  1994). 

At  night,  surface  cooling  and  associated  downward  heat  flux  cools  the  atmosphere 
leading  to  higher  atmospheric  pressure  over  the  land  surface.  The  resulting  horizontal 
pressure  gradient  drives  the  land  breeze  over  the  sea.  The  nocturnal  land  breeze  is 
thinner  and  weaker  than  the  daytime  sea  breeze  and  is  easily  prevented  by  significant 
synoptic  scale  flow.  The  land  breeze  is  rapidly  eliminated  by  convective  mixing  as  it 
flows  over  the  warmer  water  (Sun  et  al.  1998b). 

10.5.2  COASTAL  TERRAIN  SLOPES 

Most  coasts  are  characterized  by  sloping  land  surface.  These  slopes  exert  two 
independent  influences  on  local  circulations:  1)  modification  of  the  pressure  field 
through  surface  heating  on  the  sloped  terrain  and  2)  pressure  adjustments  associated 
with  topographically  forced  vertical  motions  in  a  stratified  atmosphere.  With  sunny 
conditions,  surface  heating  over  the  sloped  land  surface  causes  upslope  flow  in  the 
daytime  and  downslope  at  night  which  augments  the  sea/land  breeze  circulation  system. 
The  slope  flows  may  be  non-hydrostatic  for  steep  small  scale  slopes  and  approximately 
hydrostatic  for  larger  scale  weak  slopes  (Mahrt  1982). 

Rising  motion  forced  by  onshore  flow  and  sloped  terrain  corresponds  to  adiabatic 
cooling  in  the  stratified  atmosphere  which  in  turn  increases  the  underlying  hydrostatic 
pressure.  Much  of  this  pressure  increase  can  be  associated  with  thickening  of  a  cool 
surface  marine  layer  capped  by  an  inversion.  The  resulting  local  pressure  gradient  acts 
to  reduce  onshore  flow.  As  a  result,  onshore  flow  decelerates  before  reaching  the  coast 
and  the  local  winds  at  the  coast  are  mainly  parallel  to  the  shore  with  rising  terrain  at  the 
coast  The  parallel  flow  may  assume  the  form  of  a  low  level  jet  (Zemba,  Friehe  1987). 

Even  flow  approximately  parallel  to  the  coast  experiences  mesoscale  disturbances 
induced  by  irregularities  of  the  coastline.  Circulations  include  supercritical  flow 
(Winant  et  al.  1988;  Samelson,  Lentz  1994),  coastal-trapped  disturbances  (Holland, 
Leslie  1986;  Samelson,  Rogerson  1996;  Mass,  Albright  1987),  gravity  currents 
(Dorman  1987)  and  a  variety  of  other  mesoscale  flows  (Beardsley  et  al.  1987).  These 
mesoscale  perturbations  may  be  modulated  by  diurnal  variations  associated  with 
differential  surface  heating,  discussed  above.  In  addition,  atmospheric  and  oceanic 
circulations  are  coupled  in  the  coastal  zone.  As  one  example,  Zembe  and  Friehe  (1987) 
examine  the  influence  of  wind  acceleration  and  increased  stress,  resulting  upwelling  and 
reduction  of  surface  temperature.  The  resulting  increased  atmospheric  stability  acts  to 
reduce  the  stress. 

10.6  Conclusions 

Airflow  in  the  coastal  zone  is  complex  due  to  formation  of  internal  boundary  layers, 
diumally  varying  horizontal  pressure  gradients  and  strong  spatial  variation  of  the  wave 
field.  Assumptions  required  by  Monin-Obukhov  similarity  theory  may  not  be  met  due  to 
surface  heterogeneity,  advection  and  strong  vertical  divergence  of  the  flux.  The  surface 
airflow  is  not  only  coupled  to  the  spatially  varying  wave  field  but  also  coupled  to 
oceanic  circulations  through  the  sea  surface  temperature.  Complex  topography  at  the 
coast  further  complicate  the  total  flow  system.  The  above  interpretive  survey  did  not 
include  a  number  of  important  topics  such  coastal  zone  cloud  systems  and  special 
effects  of  bottom  bathymetry.  Improved  understanding  of  fluxes  in  the  coastal  zone 
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must  begin  with  observations  of  spatial  variation  of  the  flux  with  simple  bathymetry, 
straight  coast  line  and  no  significant  topography. 
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Summary 

Vertical  profiles  of  the  time-averaged  wind  stress,  wind  speed  and  buoyancy  flux  from  the  off-shore  tower  site 
in  the  Ris0  Air  Sea  Experiment  are  used  to  evaluate  similarity  theory  in  the  coastal  zone.  The  observed  dependence 
of  the  non-dimensional  wind  shear  on  stability  is  compared  to  the  traditional  parametrization.  Relationships 
between  the  non-dimensional  shear,  development  of  internal  boundary  layers  and  wave  state  are  explored. 

We  find  that  the  laigest-scale  turbulent  eddies  are  suppressed  in  shallow  convective  internal  boundary  layers, 
leading  to  larger  non-dimensional  shear  than  that  of  the  traditional  prediction  based  only  on  stability.  In  shallow 
stable  boundary  layers,  elevated  generation  of  turbulence  leads  to  smaller  non-dimensional  shear  compared  to 
the  traditional  prediction.  Over  young  growing  waves  in  stable  stratification,  the  observed  non-dimensional  shear 
is  less  than  that  over  older  more  mature  waves  in  otherwise  similar  conditions.  The  non-dimensional  shear  is 
a  function  of  wave  state  for  stable  conditions  even  though  the  observations  are  well  above  the  wave  boundary 
layer.  We  conclude  that  development  of  shallow  internal  boundary  layers  and  young  growing-wave  fields,  both  of 
which  are  common  in  the  coastal  zone,  can  lead  to  substantial  departures  of  die  non-dimensional  shear  from  the 
prediction  based  only  on  stability. 

Keywords:  Boundary  layer  Similarity  theory  Turbulence 


L  Introduction 


Monin-Obukhov  (MO)  similarity  theory  predicts  that  for  homogeneous  and  station¬ 
ary  conditions,  the  non-dimensional  wind  shear  (0m)  is  a  universal  function  of  atmos¬ 
pheric  stability  alone,  such  that 


0m  (z/L)  —  — — 

K*  az 


(1) 


where  z  is  height,  u  is  velocity,  z/L  is  the  stability  parameter,  L  is  the  Obukhov  length 
(Monin  and  Obukhov  1954),  which  is  a  function  of  the  friction  velocity  (k*)  and  the 
buoyancy  flux,  and  k  is  von  Karman’s  constant.  MO  similarity  theory  applies  in  the 
surface  layer  above  the  roughness  sub-layer  over  land,  and  above  the  wave  boundary 
layer  over  water  (Chalikov  and  Makin  1991).  The  top  of  the  wave  boundary  layer  has 
been  defined  as  the  height  above  the  waves  where  the  component  of  the  stress  in  phase 
with  the  waves  becomes  a  small  fraction  of  the  total  stress.  This  height  is  thought  to  be 
approximately  A/10  (Chalikov  and  Belevich  1993)  or  A/6  (Wetzel  et  al.  1996),  where  A 
is  the  wave  length  of  the  dominant  surface-wave  mode.  For  further  discussion  of  the 
wave  boundary  layer,  see  Belcher  and  Hunt  (1998),  Makin  et  al.  (1995)  and  Large 
et  al.  (1995).  The  upper  boundary  of  the  surface  layer  where  MO  similarity  theory 
applies  is  determined  by  the  assumption  that  the  height-dependence  of  the  fluxes  must  be 
neglected.  This  normally  implies  that  the  surface  layer  is  thin  compared  to  the  boundary 
layer  depth.  Grant  (1992)  has  suggested  that  MO  similarity  theory  should  be  generalized 
in  some  circumstances  to  include  the  influence  of  additional  length  scales,  such  as  the 
depth  of  the  boundary  layer. 

In  this  study,  coincident  measurements  of  vertical  profiles  of  wind  speed  and  vertical 
profiles  of  turbulent  fluxes  of  momentum  and  buoyancy,  collected  during  the  Ris0  Air 
Sea  Experiment  (RASEX),  are  used  to  evaluate  similarity  theory  hypotheses  in  the 
marine  surface  layer.  The  RASEX  dataset  is  the  first  dataset  with  detailed  vertical 
resolution  of  the  time-averaged  structure  in  off-shore  flow.  The  RASEX  off-shore  flow 

*  Corresponding  author:  College  of  Oceanic  and  Atmospheric  Sciences,  Oregon  State  University,  Corvallis,  OR 
97331,  USA. 
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observations  provide  an  opportunity  to  study  young  developing  wave  fields  and  shallow 
internal  boundary  layers  (IBLs).  Shallow  IBLs  are  common  in  the  coastal  zone  due 
to  the  roughness  and  heat  flux  discontinuity  at  the  land-sea  boundary  (Garratt  1990; 
Mahrt  et  al.  1998).  The  observed  dependence  of  the  non-dimensional  shear  on  z/L  will 
be  compared  to  the  traditional  land-based  parametrization,  and  relationships  between 
the  non-dimensional  shear,  wave  state  and  IBL  development  will  be  explored. 

2.  Similarity  theory 

The  non-dimensional  shear  of  MO  similarity  theory  has  typically  been  evaluated 
by  measuring  the  turbulent  fluxes  at  a  height  thought  to  be  in  the  surface  layer,  often 
taken  as  10  m,  and  the  shear  from  the  mean  wind  speed  at  2  levels  which  encompass 
the  turbulence  measurements.  Data  collected  over  land  in  several  field  programmes  (e.g. 
Businger  et  al.  1971;  Dyer  1974;  Hogstrom  1988)  have  been  used  to  . fit  empirically 
the  functional  form  of  <t>m(z/L).  These  studies  have  generally  confirmed  that  the 
non-dimensional  shear  appears  to  be  a  universal  function  of  z/L.  A  commonly  used 
parametrization  (Dyer  1974)  based  on  data  from  these  field  programmes  is 

<Pm(z/L )  =  {1  -  16(z/L)}-1/4;  z/L  <  0  (2) 

4>m(z/L)  =  1  +  5 (z/L);  z/L  >  0.  (3) 

An  alternative  formulation  for  the  stable  case,  proposed  by  Beljaars  and  Holtslag  (1991), 
is 

<Prn(z/L )  =  1  +  J-  \a  +  b  e~dt  (l  +  c  -  dj-)  }  (4) 

where  a  =  1, 6  =  0.667,  c  =  5,  and  d  =  0.35. 

In  local  similarity  theory,  the  turbulent  fluxes  and  stability  are  not  constrained  to 
be  constant  with  height  and  the  non-dimensional  shear  takes  the  form  (e.g.  Wyngaard 
1973;  Nieuwstadt  1984;  Sorbjan  1986) 


<Pml(z/  A)  =  —  (5) 

1  dz 

where  the  subscript  1  refers  to  the  local  value  at  a  given  height,  and  A  is  the  value  of 
L  calculated  from  the  local  height-dependent  fluxes.  This  theory  recognizes  that  local 
fluxes  should  be  related  to  the  local  shear  and  stability.  When  the  height-dependence 
of  the  fluxes  of  momentum  and  buoyancy  can  be  neglected,  local  similarity  reduces  to 
MO  similarity.  In  this  sense,  MO  similarity  is  a  subset  of  local  similarity.  Studies  that 
report  0m  based  on  measuring  fluxes  at  a  single  level,  e.g.  the  10  m  level  on  a  tower, 
are  formally  evaluating  local  similarity  theory  due  to  a  lack  of  information  on 
the  change  in  the  fluxes  with  height.  The  distinction  between  MO  and  local  similarity 
theory  becomes  uncertain  in  studies  that  do  not  evaluate  the  height-dependence  of  the 
flux.  Previous  single-level  studies  have  been  postulated  as  MO  similarity  by  implicitly 
assuming  that  the  height-dependence  of  the  flux  can  be  neglected.  This  is  probably  a 
good  assumption  in  the  lowest  few  tens  of  metres  in  the  daytime  convective  bound¬ 
ary  layer.  However,  with  internal  boundary  layers  or  morning  and  evening  boundary 
layer  transition  periods,  which  are  the  main  source  of  near-neutral  data  in  many  field 
programmes,  this  assumption  breaks  down.  Therefore,  some  observational  studies  have 
computed  stability  functions  intended  for  MO  similarity  when  in  actuality  the  stability 
functions  are  more  appropriate  to  local  similarity  theory. 
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According  to  traditional  similarity  theory,  the  non-dimensional  shear  does  not  de¬ 
pend  on  any  characteristics  of  the  surface.  Historically,  the  surface  has  been  represented 
by  the  roughness  length  (z0),  which  is  the  constant  of  integration  obtained  when  in¬ 
tegrating  Eq.  (1)  with  respect  to  z  to  obtain  (z/L)  (e.g.  Paulson  1970).  Over  solid 
surfaces,  the  roughness  length  is  usually  considered  to  be  a  constant  for  a  given  location, 
and  dependent  on  wind  direction  over  heterogeneous  surfaces.  The  roughness  length 
for  momentum  is  found  to  be  directly  related  to  the  physical  size  and  spacing  of  the 
roughness  elements  of  the  surface. 

Unlike  conditions  over  land,  the  sea  surface  roughness  field  consists  of  a  broad  band 
of  waves  which  are  generated  by  and  interact  with  the  mean  wind.  With  accelerating 
winds,  surface  gravity  waves  grow  by  extracting  momentum  from  the  mean  flow,  while 
with  decelerating  winds  and  mature  waves,  the  surface  wave  field  extracts  considerably 
less  energy  or  even  loses  energy  to  the  mean  flow.  The  wind  and  waves  in  the  coastal 
zone  are  often  in  non-equilibrium  due  to  bathometry  and  limited  fetch.  Sea  surface 
roughness  lengths  are  much  smaller  than  typical  land  values.  Over  the  sea,  the  roughness 
length  has  been  parametrized  using  the  Chamock  (1955)  relationship,  which  predicts 
that  the  roughness  length  is  proportional  to  the  wind  stress. 

Numerous  studies  have  found  that  for  a  given  wind  speed  and  stability  parameter 
(z/L),  the  drag  coefficient  and  roughness  length  are  larger  over  young  and  developing 
wave  fields  compared  with  older  wave  fields  that  are  more  in  quasi-equilibrium  with 
the  wind  (e.g.  Kitaigorodskii  1970;  Snyder  et  al.  1981;  Geemaert  et  al.  1987;  Donelan, 
1990;  Smith  et  al.  1992;  Donelan  ef  al  1993;  Vickers  and  Mahrt  1997a).  Wave  age  is 
often  used  as  a  measure  of  the  maturity  of  the  wave  field,  and  is  typically  calculated  as 
the  ratio  of  the  wave  phase  speed  of  the  dominant  wave  mode  (Cp)  (assuming  a  dominant 
mode  is  well  defined)  to  either  the  friction  velocity  or  the  wind  speed  at  some  level  near 
the  surface  (e.g.  10  m).  The  effect  of  wave  age  has  been  parametrized  within  the  existing 
framework  of  similarity  theory  by  formulating  the  roughness  length,  or  equivalently  the 
neutral  drag  coefficient,  as  a  function  of  wave  age  (e.g.  Toba  and  Koga  1986;  Maat  et  al. 
1991;  Donelan  1990;  Smith  et  al.  1992). 

Some  studies  have  approached  the  influence  of  wave  state  by  directly  examining 
the  dependence  of  <f>m  (or  equivalently  and  a  specified  z0)  on  wave  state  rather  than 
adjusting  the  roughness  length.  Bergstrom  and  Smedman  (1995)  found  that  the  non- 
dimensional  shear  was  less  with  young  waves  (small  wave-age  parameter)  compared  to 
old  waves,  but  concluded  that  the  relationship  was  not  statistically  significant,  perhaps 
due  to  their  small  range  of  conditions.  However,  their  result  is  consistent  with  enhanced 
mixing  and  stress  over  young,  growing  waves.  Davidson  (1974)  found  that  isolating 
the  dependence  of  the  non-dimensional  shear  on  wave  state  is  sensitive  to  how  the 
stability  dependence  is  removed.  He  concluded  that  for  his  open-ocean  dataset,  the  non- 
dimensional  shear  increased  with  wave  age,  in  contrast  to  traditional  MO  similarity 
theory  which  contains  no  wave  state  dependence.  In  his  dataset,  the  wave  field  was 
characterized  as  surge  swell  (values  of  Cp/u*  >  25)  representative  of  an  older  decaying 
wave  field,  conditions  were  always  stable  (z/ L  >  0),  and  almost  all  results  were  based 
on  a  single  vertical  level  (8  m).  Consistent  with  enhanced  mixing  over  younger  waves, 
Davidson  (1974)  suggested  a  form  ^m(z/L,  Cp/u»),  where  fm  decreased  with  wave 

age  for  Cp/«*  >25.  ,  , 

Miller  et  al.  (1997)  reported  observational  evidence  in  near-neutral  conditions  tor  an 

intermediate  layer  of  order  10  m  thickness  located  above  the  wave  boundary  layer  (order 
1  m)  and  below  an  outer  logarithmic  layer.  In  the  intermediate  layer,  the  wind  shear  was 
substantially  less  than  predicted  by  MO  similarity  (smaller  <pm  than  predicted).  They 
found  that  the  wave-induced  stress  associated  with  atmospheric  motions  in  phase  with 
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the  surface  waves  was  significant  in  the  underlying  wave  boundary  layer  and  was  small 
m  the  intermediate  layer. 

In  the  coastal  waters  of  the  Baltic  Sea,  Smedman  et  al.  (1995)  found  that  in  stable 
flow,  the  non-dimensional  wind  shear  and  temperature  gradient  based  on  observations 

Sea  surface  dld  not  aIways  agree  with  traditional  MO  similarity.  They 
attributed  this  to  an  extremely  shallow  surface  layer  due  to  the  presence  of  a  low-level 

Jet  P°ct0  PO  m)  and  suggested  that  the  height  of  the  jet  may  be  an  important  length 
scale.  Smedman  et  al.  (1995)  proposed  that  the  low-level  jet  suppresses  low-frequency 
motions  thus  reducing  the  transport  and  increasing  the  non-dimensional  wind  shear. 
They  found  general  agreement  with  traditional  MO  similarity  at  the  same  site  during 
similar  stability  when  no  low-level  jet  was  detected. 

Grant  (1992)  suggests  that  the  non-dimensional  gradient  0  (z/L)  for  the  near-neutral 
boundary  layer  should  be  generalized  to  be  of  the  form  <P(z/h,  h/L ,  u*/fh)  in  order 
to  indude  the  influence  of  boundary  layer  depth  (h)  on  the  non-dimensional  shear. 
Khanna  and  Brasseur  (1997)  consider  the  form  0fe/L,  h/L)  where  their  0  is  the  non- 
dimensional  gradient  of  an  arbitrary  variable.  In  the  large  eddy  simulation  results  of 
Khanna  and  Brasseur  (1997),  the  non-dimensional  shear  decreases  above  the  surface 
layer,  this  is  also  observed  by  Smedman  and  Johansson  (1997)  in  shallow  off-shore 
boundary  layers.  Apparently,  the  vertical  gradient  decreases  more  rapidly  with  height 

*  wi,?  d0es  flux*  No  spedfic  formulation  for  a  generalization  of  0 

to  include  the  dependence  on  boundary  layer  depth  has  been  attempted.  If  the  influence 

o  oundary  layer  depth  on  0m  extends  downward  to  the  surface,  or  if  the  influence  of 
wave  state  on  0^  extends  upward  above  the  wave  boundary  layer,  then  both  local  and 
similarity  theory  break  down.  We  intend  to  consider  several  possibilities:  (a)  MO 
W1?  extstog  stability  functions  is  valid;  (b)  MO  similarity  theory  is  valid  but 
the  traditional  stability  functions  require  modification;  (c)  MO  similarity  is  not  valid  but 

' m?  1S/??J?PnT  Wlth  possible  modification  of  the  stability  functions;  and 
(d)  both  local  and  MO  similanty  theory  are  invalid  and  additional  information  such  as 
boundary  layer  depth  and  wave  state  are  required. 


3.  Dataset 

The  full  RASEX  instrumentation  is  described  in  Barthelmie  et  al  (1994)  and 
Hdjstrup  et  ^(1997^  Data  used  in  this  study  were  collected  at  the  sea  mast  west 
tower  located  2  km  off  the  coast  of  Lolland,  Denmark  in  4  m  of  water.  Observations 
taken  from  April  through  November  1994,  including  the  intensive  campaigns  during  the 
spring  and  autumn,  are  considered.  The  off-shore  tower  site  is  representative  of  an  inland 
sea  or  large  lake,  with  only  small  amplitude  swell  and  predominantly  wind  driven  waves 
ampIltudes  ,ess  than  1  m  and  wavelengths  less  than  16  m  (Vickers 

i  9  h  1  f°nS  °n'Sh0re  floW’  ^  wave  boundary  layer  may  sometimes 
include  the  lowest  turbulent  flux  measurement  level  at  3  m  (Hare  et  al.  1997)  The  lowest 

observational  level  used  in  this  study  (6  m)  is  expected  always  to  be  above  the  top  of  the 
wave  boundaty  layer  based  on  the  maximum  observed  wavelength.  Shallow  IBLs  are 
common  in  off-shore  flow  where  the  fetch  ranges  from  2  to  5  km  (Mahrt  et  al.  1998) 

On-shore  flow,  where  the  fetch  ranges  from  15  to  25  km,  is  still  fetch-limited  compared 
to  the  open  ocean.  F 

Fast-response  (20  Hz)  observations  from  Gill/Solent  sonic  anemometers  (hereafter 

^hShPi°Sltl°ne^  6’  10’  18  and  32  m  above  the  sea  surface  were  used  to  estimate 
the  turbulent  wind  stress  and  buoyancy  flux  profiles.  Turbulence  measurements  from 
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the  instrument  at  the  3  m  level  were  excluded  in  this  study  due  to  potential  flow- 
distortion  problems.  Fluxes  calculated  from  the  45  m  level  sonic  data  were  found  to  be 
erratic  and  unreliable,  presumably  due  to  an  instrumental  problem,  and  were  excluded. 
The  tilt-correction  method  used  for  the  sonics  was  presented  in  Mahrt  et  al.  (1996). 
The  buoyancy  flux  was  calculated  from  the  sonic  virtual-temperature  fluctuations,  after 
applying  a  correction  for  acoustic  wave  bending  that  can  be  important  in  strong  winds. 
Vertical  profiles  of  wind  speed  were  measured  using  cup  anemometers  (P224b  sensor) 
located  7,  15,  20,  29  and  38  m  above  the  sea  surface.  Time  periods  where  the  cup 
anemometer  wind  speed  was  less  than  2.5  m  s-1  were  excluded;  because  the  cup 
anemometers  were  calibrated  for  the  generally  strong  winds  at  the  site,  readings  of 
less  than  2.5  m  s-1  are  not  reliable.  Temperature  profiles  were  measured  using  a 
thermometer  at  10  m  (Risp  PT  100),  and  temperature-difference  sensors  (Ris0  PT  500) 
to  provide  the  vertical  gradient  of  temperature  over  the  layers  10  to  24  m  and  10  to 
47  m.  The  fluctuating  wave  height  was  measured  using  an  acoustic  wave  recorder.  A 
characteristic  wave  amplitude,  wave  phase  speed  and  wavelength  were  calculated  as 
described  in  Vickers  and  Mahrt  (1997a).  All  time  periods  with  wind  direction  from  the 
sector  including  339  deg.  through  106  deg.  were  excluded,  because  of  potential  tower 
interference  and  possible  flow  distortion  due  to  wind  turbines  located  north-east  of  the 

tower.  , .  .  , 

The  sonic  and  cup  anemometer  and  wave  recorder  measurements  were  subjected 

to  quality  control  screening  as  described  in  Vickers  and  Mahrt  (1997b).  The  automated 
quality  control  procedure  identifies  periods  with  potential  problems  with  instruments  or 
data  recording.  The  package  of  tests  checks  for  excessive  electronic  spiking,  inadequate 
amplitude  resolution,  signal  dropout,  unrealistic  magnitude,  extreme  higher  moment 
statistics  and  near  discontinuities  in  the  first  two  moments.  Visual  inspection  of  all  data 
flagged  by  the  automated  procedures  was  carried  out,  either  to  verify  an  instrument 
problem  or  to  identify  unusual  but  plausible  atmospheric  behaviour.  In  the  former  case, 
the  flag  was  verified  and  the  data  were  eliminated. 

Direct  observations  of  IBL  depth  at  the  site  were  not  routinely  available  during  RA- 
SEX.  Profiles  of  potential  temperature,  humidity  and  winds,  measured  by  tethersondes 
released  from  land  in  the  vicinity  of  the  site  during  12  different  days  in  October  and 
early  November  1994,  were  examined.  A  parametrized  IBL  depth  (e.g.  H0jstrup,  1981) 

h  =  Caw(X/u ),  (6) 

where  ow  is  the  vertical  velocity  fluctuation,  X  is  the  fetch  distance  and  C  is  an  empirical 
coefficient,  was  found  to  correlate  reasonably  well  with  the  IBL  depth  as  inferred  from 
the  tethered  balloon  soundings  for  the  limited  number  of  comparisons  available.  We 
found  that  setting  C  =  0.4  gave  a  better  fit  to  the  data  than  the  original  coefficient  of 
unity,  and  use  C  =  0.4  in  all  further  calculations.  Using  the  vertical  velocity  fluctuations 
from  the  18  m  sonic,  the  38  m  wind  speed  from  the  cup  anemometer  and  the  actual 
fetch  distance,  the  parametrized  h  at  the  tower  location  varies  from  a  minimum  of  20  m 
in  off-shore  flow  (1:100  stable  IBL  growth  rate)  to  a  maximum  of  400  m  in  on-shore 
flow  (1:30  unstable  IBL  growth  rate).  We  note  that  u  in  Eq.  (6)  should  represent  the 
mean  wind  at  the  top  of  the  IBL,  which  results  in  an  iterative  expression  for  h.  Here,  we 
approximate  u  at  height  h  using  the  38  m  wind  speed. 

(a)  Fluxes 

Eddy  correlation  fluxes  of  momentum  and  virtual  temperature  were  calculated  from 
sonic  anemometers  at  4  levels  (6,  10,  18  and  32  m)  to  estimate  profiles  of  the  wind 
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stress  and  Obukhov  length.  An  averaging  time-scale  of  10  minutes  was  selected  to  define 
e  turbulent  fluctuations,  and  the  flux  computed  from  such  fluctuations  was  averaged 
over  one-hour  periods.  Therefore,  the  fluxes  include  all  motions  on  time-scales  of  10 
£r  “n  t0  the  smallest  resolvable  scales  determined  by  the  sampling 
rate  (20  Hz).  These  averaging  times  were  chosen  to  optimize  the  correlation  between 
the  mean  wind  speed  and  the  wind  stress  (Mahrt  et  al.  1996).  Averaging  times  longer 
than  one  hour  incorporate  too  much  non-stationarity,  and  shorter  averaging  times  do 
not  include  enough  samples.  Elimination  of  records  with  large  non-stationarity  and  flux 
sampling  errors  are  discussed  in  the  appendix.  The  wind  stress  was  calculated  as  the 
vector  magnitude  of  the  one-hour  averages  of  the  stress  components. 

(b)  Wind  profiles 

The  wind  shear  was  estimated  by  applying  the  method  of  least  squares  to  fit  a 
second-order  polynomial  in  In (z)  to  the  profiles  of  one-hour  mean  wind  speed  from 
die  cup  anemometers  at  5  levels  (7,  15,  20,  29  and  38  m),  and  also  a  derived  level  at 
toe  roughness  height  (za)  where  the  wmd  speed  was  assumed  to  be  zero.  The  precise 
heights  of  the  wmd  speed  measurements  were  adjusted  to  account  for  the  fluctuating 
mean  sea  level.  The  wmd  speeds  from  the  cup  anemometers  were  re-calibrated  using 
the  post-expenment  calibration  coefficients  (Hpjstrup,  personal  communication).  The 
wind  shear  was  evaluated  analytically  at  the  flux  measurement  heights  (6,  10,  18  and 

„  second-°rder  polynomial  fit  to  the  wind  speed  profile 

^  Hogstron\(1988).  The  Chamock  relationship  with  a  coefficient  of  0.017  (Garratt 
1977)  was  used  as  a  lower-boundaiy  condition  for  the  fit.  The  fit  to  the  wind  speed 
profile  was  not  sensitive  to  the  choice  of  this  coefficient  nor  application  of  a  wave-age 
dependent  Chamock  relationship.  6 

•  3  STSi?Vity  s?dy’  *e  wind  sPeed  Profiles  were  fitted  using  the  cup  anemometer 

wmd  spee<is  alone  the  sonic  speeds  alone  and  a  combined  set  of  cup  and  sonic  wind 
speeds  We  also  calculated  the  shear  using  finite  differences  applied  to  adjacent  cup 
levels,  differencing  with  respect  to  both  linear  z  and  ln(z).  In  most  instances,  the  profile 
of  wind  speed  from  the  cups  was  smoother  than  both  the  profile  from  the  sonics  and 
nt  Secfr01?  the  <;ombmed  cuP  and  sonic  set.  In  terms  of  variance  explained,  the 
Pmd  p™filesuwere  better  represented  by  the  second-order  polynomial  than 
were  die  profiles  from  the  some  speed  measurements  or  the  combined  cup  and  sonic 
set.  The  root-mean-square  error  of  the  model  fit  to  the  wind  speed  profile  evaluated 

at  the  5  cup  anemometer  levels  ranged  from  0.03  m  s"1  to  0.19  m  s'*  and  averaged 
0.09  m  s  .  6 


4.  Non-dimensional  shear 

The  observed  dependence  of  0mi  on  z/A  for  all  the  data  including  all  4  vertical  levels 
and  both  on-shore  and  off-shore  flow  regimes  is  shown  in  Fig.  1.  The  observed  <bmX  is 
less  than  the  traditional  prediction  (Eqs.  (2)-(4))  for  the  entire  range  of  stability  (greater 
stress  re^tivetothe  shear).  The  near-neutral  (-0.01  <  z/A  <0.01)  value  of/mlis  0.95 

^  Standard  erTOr  of  *e  avera§e  0ml  for  z/A  categories  is  relatively 
small  for  the  near-neutral  categories,  and  larger  for  unstable  and  stable  categories  where 
the  number  of  observations  is  less. 

For  estimating  0m  for  MO  similarity  at  the  10  m  level,  we  excluded  cases  with 
die  largest  vertical  flux  divergence  based  on  the  fluxes  at  the  6,  10  and  18  m  levels  This 
flux  divergence  criterion  excluded  approximately  one-half  of  the  data.  For  the  remaining 
cases  used  m  calculating  <pm,  the  magnitude  of  the  stress  divergence  for  the  layer  from 
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Figiire  1.  The  mean  (circles)  and  plus  and  minus  one  standard  error  of  the  local  non-dimensional  shear  W>mi)  for 

stability  (z/A)  categories  (see  text)  for  all  the  data.  The  solid  curve  is  a  fit  to  the  data:  <pm\  =  U  ”  ’ 

z/A  <  0,  and  0ml  =  1+  2.5(z/A);  z/A  >  0,  where  A  is  the  Obukhov  length  based  on  local  fluxes  The  jiash- 
dot  curve  is  the  traditional  Dyer  (1974)  parametrization.  The  dashed  curve  is  the  Beljaars  and  Hoitslag  (1991) 

formula. 


6  m  to  18  m  was  always  less  than  8  x  10“3  m  s-2  and  averaged  5x10  ms  .As 
we  found  for  local  similarity,  the  non-dimensional  shear  of  MO  similarity  theory  (<pm)  is 
also  less  than  the  traditional  prediction  (Fig.  2).  The  mean  <pm  at  10  m  for  weakly  stable 
conditions  is  larger  than  <p mi,  closer  to  both  the  predictions  (Eqs.  (3H4))  and  srmilar  to 
the  estimate  of  0ml  when  excluding  off-shore  flow  cases  (Fig.  3).  This  reflects  the  fact 
that  large  flux  divergence  is  common  in  off-shore  stable  flow.  In  this  regard,  eliminating 
periods  with  large  flux  divergence  is  similar  to  eliminating  periods  with  stable  off-shore 
flow.  In  unstable  conditions,  the  observed  10  m  <f)m  is  less  than  4>m\,  further  from  the 
prediction  (Eq.  (2))  and  similar  to  <t>m i  when  excluding  shallow  IBL  cases  (Fig.  4).  For 
the  unstable  case,  eliminating  periods  with  large  flux  divergence  also  tends  to  eliminate 
periods  with  shallow  IBL  depth.  We  conclude  that  eliminating  convective  cases  with 
large  flux  divergence  does  not  improve  the  agreement  between  the  observed  (f>m  and  the 

traditional  prediction  (Eq.  (2)).  . 

Note  that  the  range  of  the  stability  parameter  for  the  <pm  dataset,  where  z 10  m,  is 
considerably  less  than  for  the  <f>mi  dataset  which  includes  z  =  6,  10,  18  and  32  m.  The 
standard  error  bars  are  larger  for  than  for  <pm\  partly  due  to  the  smaller  sample  size. 

The  mean  $mi(z/A)  partitioned  into  on-shore  and  off-shore  flow  categories  is  shown 
in  Fig.  3.  In  this  dataset,  the  majority  of  the  strongly  convective  cases  are  found  in 
on-shore  flow  while  the  majority  of  the  strongly  stable  cases  are  found  in  off-shore 
flow.  This  complicates  attempts  to  isolate  the  effects  of  stability  and  fetch.  The  range 
of  z/A  for  an  individual  stability  category  was  expanded  in  Fig.  3  compared  to  Fig.  1, 
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Figure  2.  The  mean  and  standard  error  of  the  non-dimensional  shear  (tf^)  for  10  m  data  only  for  stability  {z/L) 
categories  (see  text).  The  solid  curve  is  a  fit  to  the  local  non-dimensional  shear  (p^  (as  in  Fig.  1).  The  dash-dot 
curve  is  the  traditional  Dyer  (1974)  parametrization.  The  dashed  curve  is  the  Beljaars  and  Holtslag  (1991)  formula. 


to  ensure  an  adequate  number  of  samples  in  each  stability  category  when  partitioning 
the  data  by  both  stability  and  flow  regime.  For  off-shore  neutral  and  stable  flow,  <f> mi 
is  significantly  less  than  predicted  by  previous  relationships.  In  the  transition  region 
downstream  from  a  rough-to-smooth  roughness  change,  the  mean  shear  may  adjust  more 
rapidly  to  the  smoother  ocean  surface  compared  to  the  stress,  preventing  equilibrium 
between  the  stress  and  the  shear  (Peterson  1969).  As  a  result,  the  non-dimensional  shear 
may  be  smaller  in  the  coastal  zone  than  that  predicted  by  traditional  similarity  theory. 
A  mechanism  to  account  for  the  difference  between  the  on-shore  and  off-shore  flow 
unstable  cases  is  discussed  next. 


(a)  IBL  depth 

The  stability  dependence  of  </>mi  for  two  z/h  categories  is  shown  in  Fig.  4  where  h 
is  the  IBL  depth  estimated  from  Eq.  (6).  The  largest  values  of  z/h  (small  h)  occur  in 
off-shore  flow  due  to  the  shorter  fetch,  and  partitioning  the  data  by  z/h  also  tends  to 
partition  the  data  into  off-shore  and  on-shore  flow  regimes.  Therefore,  other  factors  in 
addition  to  IBL  depth  may  cause  differences  between  the  two  z/h  classes. 

In  convective  IBLs,  <pm\  increases  with  increasing  z/h  with  the  exception  of  the 
most  unstable  category,  that  also  has  the  least  amount  of  data.  We  attribute  the  increase 
of  0mi  with  increasing  z/ h  to  suppression  of  the  largest  convective  eddies  by  the  shallow 
IBL  top  (Mahrt  et  al.  1998),  which  acts  to  reduce  the  stress  relative  to  the  shear.  Multi¬ 
resolution  spectra  (Howell  and  Mahrt  1997)  of  the  wind  stress  from  the  6  and  32  m  levels 
support  this  theory  in  that  lower-frequency  (large  eddy)  transport  is  indeed  suppressed 
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Figure  3.  The  mean  local  non-dimensional  shear  (0mi)  for  stability  (z/ A)  categories  (see  text)  for  on-shore  flow 
(solid)  and  off-shore  flow  (dotted).  The  dash-dot  curve  is  the  traditional  Dyer  (1974)  parametrization.  The  dashed 

curve  is  the  Beljaars  and  Holtslag  (1991)  formula. 


by  shallow  EBL  depth  (Fig.  5).  The  multi-resolution  spectra  orthogonally  decompose  the 
flow  into  simple  unweighted  averages  for  different  averaging  windows.  These  spectra 
yield  results  similar  to  Fourier  spectra  but  satisfy  Reynolds  averaging  at  all  scales  and 
do  not  assume  periodicity. 

In  contrast  to  unstable  conditions,  the  spectra  for  neutral  and  stable  conditions 
show  litde  difference  in  the  scale-dependence  of  the  stress  over  a  wide  range  of  z/h. 
Our  hypothesis  is  that  in  neutral  and  stable  conditions,  the  eddies  are  too  small  to  be 
suppressed  by  the  shallow  EBL  depth. 

The  agreement  with  traditional  similarity  for  convective  EBL  cases  is  better  when 
combining  on-shore  and  off-shore  flow  (Fig.  1),  than  when  partitioning  the  data  by  either 
flow  regime  or  z/h  (Figs.  3  and  4).  This  agreement  is  somewhat  circumstantial.  In  on¬ 
shore  flow,  where  the  IBL  depth  is  typically  large  compared  to  the  observation  height, 
and  where  we  would  expect  the  best  agreement  with  the  traditional  parametrization,  <j>m\ 
is  less  than  predicted.  In  off-shore  flow  with  shallow  convective  IBLs,  <f>m \  is  larger  than 
predicted.  These  data  suggest  that  for  the  convective  case,  <pmi  is  less  than  predicted  by 
traditional  similarity  but  increases  with  increasing  z/h. 

In  stably  stratified  IBLs,  0mi  decreases  with  increasing  z/h  (Fig.  4).  That  is,  shallow 
stable  IBLs  lead  to  relatively  large  stress  and  therefore  smaller  <pm\  than  predicted.  We 
attribute  this  to  the  top-down  nature  of  the  turbulence  in  the  presence  of  shallow  stable 
IBLs  in  advective  conditions  where  the  stress  decreases  with  height  more  slowly  than 
the  mean  shear,  or  where  it  may  even  increase  with  height  (Fig.  6).  This  is  apparently 
due  to  suppression  of  turbulence  near  the  surface  by  the  stratification  in  concert  with 
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the  parametrized  IBL.  The  dash-dot  curve  is  the  traditional  Dyer  (1974)  parametrization.  The  dashed  curve  is 

the  Beljaars  and  Holtslag  (1991)  formula. 


shear-generation  of  turbulence  aloft,  although  horizontal  advection  of  turbulence  may 
ahso  be  important.  The  cornposited  wind  profiles  for  this  case  (not  shown)  correspond 
to  increasing  mean  shear  with  height.  The  elevated  shear  may  be  associated  with  flow 
accderation  resulting  from  partial  decoupling  from  the  sea  surface,  sometimes  observed 
with  the  flow  of  warmer  air  over  cooler  water  (Smedman  et  al.  1995).  When  the  stress 
ncreases  with  height  the  usual  atmospheric  boundaiy  layer  is  turned  upside  down,  in 
the  primary  source  of  the  turbulence  is  aloft  rather  than  at  the  surface.  The  physical 
interpretation  of  the  IBL  depth  is  ambiguous  in  such  cases. 

fnr  JHmperatTe  advection  is  also  important  in  off-shore' flow.  Non-zero  heat  flux  is 
forced  pnmanly  by  temperature  advection.  That  is,  the  heat  flux  divergence  is  balanced 
ofnJl  P/  5a^re  adve<m<>n  since  the  local  time-change  of  temperature  is  usually  small. 

Thertw  depe?dK  T°n  ?e  hef  flux’  Z/L  is  correlated  with  temperature  advection. 
rsiiv  ™  ?’  any  stability-dependent  non-dimensional  quantities  over  the  sea  automati¬ 
cally  contain  advective  effects,  particularly  in  the  coastal  zone. 


(b)  Wave  age 

Fie  T7  CYonnolty  dePendence  of  fml  for  two  wave  age  (Cp/m*)  regimes  is  shown  in 
tig.  7.  Young  growing-wave  fields  preferentially  occur  with  off-shore  flow,  but  are 

!  !n  fi°m^0n  yuf0Und  10  Stf0ng  on*shore  flow.  In  stable  conditions,  0ml  increases 

agC  a$  aIS°  °bSerVed  by  Bergstr6m  aad  Smedman 
,1995)  and  Davidson  (1974).  The  self-correlation  between  wave  age  and  <pm ]  (both 
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Figure  5.  Average  multi-resolution  spectra  of  the  one-hour  average  kinematic  wind  stress  (m  s  x  100)  at  the 
6  m  level  in  unstable  conditions  for  three  z/h  categories  (see  text):  (i)  z/ ft  =  0.1  (solid);  (ii)  z/h  -  0.02  (dashed); 
and  (iii)  z/h  =  0.01  (dotted)  where  h  is  the  parametrized  internal  boundary  layer  depth.  The  time-scale  axis  is 

analogous  to  the  period  in  Fourier  analysis. 


contain  «*)  may  influence  the  result.  The  self-correlation  between  z/L  and  <pm  can  also 
be  large  (L  contains  u*).  We  have  attempted  to  reduce  the  self-correlation  by  jointly 
stratifying  <pm j  by  both  z/A  and  Cp/u*  in  Fig.  7.  In  a  separate  calculation,  we  partitioned 
the  data  into  Cp/ut  wave-age  regimes,  where  this  wave-age  parameter  is  based  on  ut, 
the  estimated  wind  speed  at  the  top  of  the  wave  boundary  layer,  instead  of  «*,  and  thus 
is  not  subject  to  self-correlation.  This  calculation  confirmed  a  real  dependence  of  <pm\  on 
wave  state.  We  conclude  that  in  stable  conditions,  the  mean  shear  does  not  completely 
adjust  to  the  growing-wave  field  leading  to  a  dependence  of  (pm\  on  wave  age. 

The  possible  dependence  of  <pmi  on  wave  age  contradicts  the  assumption  of  simi¬ 
larity  theory  that  the  local  flux-gradient  relationship  in  the  surface  layer  is  independent 
of  surface  conditions.  In  traditional  similarity  theory  for  the  surface  layer  (above  the 
wave  boundary  layer),  the  influence  of  surface  characteristics  is  included  in  the  rough¬ 
ness  length  which  appears  in  the  vertical  integration  of  4>m\  downward  to  the  surface.  In 
this  dataset,  the  observation  heights  are  well  above  the  estimated  wave  boundary  layer 
height,  so  that  dependence  of  (pm\  on  wave  state  is  not  expected  on  the  basis  of  traditional 
similarity  theory. 

For  unstable  conditions,  0ml  is  observed  to  be  larger  over  younger  waves,  contrary 
to  the  result  found  for  stable  flow  (Fig.  7).  This  may  be  due  to  the  dominant  influence 
of  shallow  IBL  depth  rather  than  wave  state  in  modifying  <f>m\  for  convective  cases 
(section  4(a)).  It  may  also  be  related  to  correlation  between  wave  age  and  other  processes 
related  to  flow  direction,  since  the  difference  between  <pm\  for  younger  waves  compared 
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Figure  6.  Average  profiles  of  the  kinematic  wind  stress  (m2s2  x  100)  and  drag  coefficient  (Crf  x  1000)  in  stable 

SST  r5, '  \°-3).SL*£ 200  m  »<!  *  <  200  m  (dasL)  where™  i'  VSlrnS 

internal  boundary  layer  depth  and  Obukhov  length  L  is  calculated  using  the  average  fluxes  from  the  6  and  10  m 

levels.  The  average  stability  (z/ L)  and  wind  speed  profile  for  both  h  categories  are  nearly  identical. 


with  older 
(Fig.  3). 


waves  is  similar  to  the  difference  when  partitioning  the  data  by  flow  direction 


5.  New  formulation 

With  off-shore  flow  in  the  coastal  zone,  shallow  IBL’s  are  common  due  to  short 
fetch.  Significant  vertical  stress  and  buoyancy  flux  divergences  are  typical  due  to  shallow 
IBL  depth  and  strong  temperature  advection.  The  wind  and  waves  can  be  in  non- 
equilibnum  due  to  the  bathometry  and  limited  fetch.  Relatively  large  stress  generated 
over  land  may  be  advected  over  the  water.  The  mean  flow  normally  accelerates  due  to 
smaller  surface  roughness  over  the  water.  All  of  these  mechanisms  influence  the  non- 
dimensional  shear. 

We  expect  that  the  conditions  required  for  application  of  similarity  theory  will  be 
best  satisfied  when  advective  effects  are  minimal  (on-shore  flow)  and  the  observation 
height  is  sufficiently  below  the  IBL  depth  (z/h  <  1/4).  The  stability  dependence  for 
data  satisfying  these  conditions  is  shown  in  Fig.  8  and  a  fit  to  the  data  is: 

0ml  =  {1  -  35(z/A)}“1/4;  z/A  <  0  (7) 

0ml  =  {1  +  16(z/A)}1/3;  z/A  >  0.  (8) 

For  this  reduced  dataset,  a  von  Karman  constant  (k)  equal  to  0.39  is  required  to  make 
0ml (z/ A  =  0)  equal  to  unity,  and  the  coefficients  in  Eqs.  (7)-(8)  are  based  on  k  =  0.39. 
The  model  predicts  unstable  values  of  0ml  that  are  approximately  15-20%  smaller 
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Figure  7.  The  mean  local  non-dimensional  shear  (0m i)  for  stability  (z/A)  categories  (see  text)  for  two  wave- 
age  (Cp/tt.)  categories:  (i)  Cp/u,  <  15  (solid);  and  (ii)  Cp/«„  >  15  (dotted),  where  Cp  is  the  phase  speed  of  the 
dominant  wave  mode.  The  dash-dot  curve  is  the  traditional  Dyer  (1974)  parametrization.  The  dashed  curve  is  the 

Beljaars  and  Holtslag  (1991)  formula. 


than  the  Dyer  (1974)  formula  (greater  stress  relative  to  the  shear).  It  is  not  clear  if  the 
relatively  weak  shear  observed  here  is  related  to  the  weak  shear  layer  observed  by  Miller 
et  al.  (1997).  For  stable  conditions,  Eq.  (8)  predicts  values  of  <pm\  significantly  less  than 
Dyer  (1974)  for  values  of  z/A  exceeding  0.2.  The  deviation  of  the  non-dimensional 
shear  from  Dyer’s  formula  is  in  the  same  direction  as  the  Beljaars  and  Holtslag  (1991) 
correction,  but  much  greater  in  magnitude.  Howell  and  Sun  (1999)  further  document  the 
overprediqtion  of  <pm  by  traditional  stability  functions  for  strong  stability,  although  the 
deviations  from  Dyer  are  greater  in  the  present  study. 

To  better  isolate  the  potential  influence  of  IBL  depth  and  wave  age  on  <pm\,  we 
calculated  the  residuals  from  the  model,  defined  as  the  observed  (pm\  minus  Eqs.  (7)-(8), 
and  partitioned  them  by  z/h  and  wave-age  categories  for  unstable  and  stable  periods. 
The  Observed  <pm\  dataset  includes  all  the  on-shore  and  off-shore  flow  periods. 

For  unstable  conditions,  the  <j>m\  residuals  increase  with  increasing  z/h  for  z/h  < 
1/3  for  all  wave-age  regimes  except  the  youngest  wave-age  category  (Table  1).  The 
main  effect  here  appears  to  be  due  to  suppression  of  large  convective  eddies  by  shallow 
IBL  depth  (section  4(a)).  For  very  shallow  IBL  depth  {z/h  >  1/3),  (pm\  decreases  with 
z/h.  We  attribute  this  to  the  elevated  generation  of  turbulence.  In  unstable  flow,  there  is 
no  clear  relationship  between  the  <pm\  residuals  and  wave  age. 

The  dependence  of  the  unstable  <pm\  residuals  on  z/h  suggests  including  the  depen¬ 
dence  of  the  non-dimensional  shear  on  z/h  for  convective  IBLs.  Here  we  choose  the 
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(1 974)  pSie2itior?he C°RST  ^=°-39' 11,6  curve  *  *e  traditional  Dyer 

(  )  Parametnzation.  The  dashed  curve  is  the  Beljaars  and  Holtslag  (1991)  formula. 


L  ^'^nsional  W,ND  SHEAR  Residuals  (observed  minus 

TH? tows Irf  wl  rSn0LUT  ARE  WAVE'AGE  (Cp/«»)  CATEGORIES  AND 
HE  ROWS  ARE  z/h  CATEGORIES.  NUMBERS  OF  OBSERVATIONS  ARE  SHOWN  IN 


PARENTHESIS 


5-10 

10-15 

15-20 

20-25 

25-30 

0.32-1 

0.1-0.32 

0.032-0.1 

0.01-0.032 

-0.08  (20) 
0.01  (41) 
0.15(08) 
0.14(07) 

0.21  (13) 
0.30  (40) 
0.03  (39) 
-0.03  (29) 

-0.06  (14) 
0.00  (43) 
-0.02  (63) 
-0.07  (34) 

-0.10(09) 
-0.03  (17) 
0.01  (17) 
-0.09  (10) 

0.19(04)' 
-0.05  (06) 

TABLE  2. 

As  Table  1  but  for  z/a  >  0 

5-10 

10-15 

15-20 

20-25 

25-30 

0.316-1 

0.1-0.316 

0.0316-0.1 

0.01-0.0316 

-0.56  (43) 
-0.11  (39) 
-0.03  (28) 
-0.08(16) 

-0.29  (43) 
-0.06  (62) 
-0.02  (52) 
0.17(19) 

-0.17(26) 
-0.03  (30) 
0.13(21) 
0.37  (10) 

0.23  (08) 
0.07(18) 
0.22(21) 
0.25(10) 

0.80(10) 
1.42(06) 
0.47  (05) 
0.44  (03) 
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form  for  the  unstable  case  as 

<t>m i  =  {1  -  35(z/ A)}- 1/4  +  (~z/A)ll2(z/h)(l  -  1-5 z/h).  (9) 

The  term  (z/ft)(l  -  1.5 z/h)  augments  <pm\  from  the  stability-dependence  prediction 
for  the  parameter  range  0  <  z/h  <  2/3,  and  reduces  <pm\  for  z/h  >  2/3.  We  include  this 
expression  mainly  to  clarify  the  observed  z/h  dependence  of  <pmi  in  convective  IBLs. 
A  more  precise  calibration  of  this  expression  would  require  a  larger  dataset  and  more 
explicit  determination  of  the  boundary  layer  depth. 

For  stable  conditions,  the  0mi  residuals  clearly  increase  with  increasing  wave  age 
for  the  entire  range  of  z/h  (Table  2).  Theoretically,  the  shear  in  the  surface  layer 
should  adjust  to  the  stress  such  that  4>m\  is  independent  of  wave  age.  Apparently 
with  stable  stratification,  the  adjustment  time-scale  of  the  weak  turbulence  may  be 
too  slow  to  maintain  equilibrium  with  the  changing  surface  state.  Consequently,  the 
assumptions  required  for  similarity  theory  are  not  met.  As  a  possible  result,  the  flux- 
gradient  relationship  (<£mi)  at  the  observational  level  now  depends  on  wave  state.  A 
simple  model  which  demonstrates  the  observed  dependence  of  the  <pm \  residuals  on  wave 
age  for  the  stable  case  is 

•  0ml  =  (1  +  16(z/A)}1/3  +  (z/A)2(Cp/w*  —  18)/20,  (10) 

where  the  coefficients  (18  and  20)  were  determined  using  the  average  residuals  for  the 
wave-age  categories  in  Table  2.  The  wave-age  term  augments  the  stability-dependent 
prediction  of  4>m\  for  older-wave  fields  (Cp/m*  >  18)  and  reduces  <pm\  over  younger 
waves.  This  model  predicts  zero  wave-age  influence  for  wave  age  equal  to  18,  which  is 
smaller  than  the  value  of  25  thought  to  represent  minimal  wind-wave  coupling  (Volkov 
1970). 

In  stable  flow  over  young  waves  (Cp/w*  <  20),  the  <pm\  residuals  clearly  become 
more  negative  with  increasing  z/h  (stronger  stress  relative  to  the  shear),  while  over  older 
waves  the  dependence  on  z/h  is  weak  (Table  2).  The  decrease  in  <pm\  with  increasing 
z/h  in  stable  flow  was  noted  in  section  4(a),  and  attributed  to  a  possible  top-down  nature 
of  the  turbulent  transport.  Table  2  shows  that  the  combined  influence  of  young  waves 
and  shallow  stable  IBL  depths  leads  to  the  largest  negative  <f> mi  residuals  (smallest  non- 
dimensional  shear  compared  to  Eqs.  (7)— (8))  in  this  dataset.  In  this  situation,  the  stress 
is  large  relative  to  the  shear  due  to  the  influence  of  wind-wave  coupling  and  the  elevated 
generation  of  turbulence.- 


6.  Conclusions 

For  convective  internal  boundary  layers  in  the  coastal  zone,  the  non-dimensional 
wind  shear  is  less  than  predicted  by  traditional  similarity  but  increases  with  increasing 
z/h.  The  largest  convective  turbulent  eddies  are  suppressed  by  the  top  of  the  internal 
boundary  layer,  leading  to  an  increase  in  the  non-dimensional  shear  with  increasing 
z/h. 

In  stable  internal  boundary  layers,  elevated  generation  of  turbulence  can  lead  to 
an  upside  down  boundary  layer  where  the  stress  may  even  increase  with  height.  As  a 
result,  the  non-dimensional  shear  decreases  with  increasing  z/h.  This  is  apparently  due 
to  suppression  of  turbulence  near  the  surface  by  the  stratification  and  shear-generation 
of  turbulence  aloft.  However,  advection  of  stress  may  also  be  important. 

In  stable  stratification,  smaller  non-dimensional  shear  is  observed  over  young  grow¬ 
ing  waves  compared  to  older  waves.  We  speculate  that  with  the  long  adjustment  time- 
scale  and  weaker  turbulence  associated  with  stable  stratification,  the  mean  shear  in  the 


75 


D.  VICKERS  and  L.  MAHRT 


coastal  zone  does  not  completely  adjust  to  the  wind-wave  coupling  influence  on  the 
stress.  In  addition,  the  surface  waves  may  effectively  perturb  the  adjacent  stable  atmos¬ 
phere,  since  vertical  velocity  fluctuations  in  a  stratified  atmosphere  induce  temperature- 
related  pressure  perturbations  (such  as  internal  gravity  waves).  Therefore,  the  non- 
dimensional  shear  depends  on  wave  state,  even  though  the  observational  levels  are  well 
above  the  wave  boundary  layer  where  individual  eddies  are  coherent  with  surface  waves. 
For  convective  and  neutral  conditions,  wave  state  does  not  appear  to  have  an  influence 
on  the  non-dimensional  shear. 

A  formulation  for  the  non-dimensional  shear  (Eqs.  (7)-(8))  in  the  marine  surface 
layer  is  constructed  as  a  function  of  the  traditional  stability  parameter,  based  on  data 
where  the  conditions  required  for  application  of  similarity  theory  are  expected  to  be  met. 
Compared  to  the  traditional  land-based  parametrization,  the  new  formulation  predicts 
slightly  smaller  non-dimensional  shear  for  unstable  conditions  and  significantly  smaller 
non-dimensional  shear  for  strong  stability.  The  formulation  is  generalized  to  include 
the  influence  of  internal  boundary  layer  depth  (Eq.  (9))  and  wave  state  (Eq.  (10)). 
This  formulation  requires  evaluation  from  a  more  extensive  dataset  to  better  isolate  the 
influences  of  stability,  internal  boundary  layer  depth  and  wave  state,  that  are  sometimes 
inter-correlated.  The  formulation  of  the  non-dimensional  shear  has  been  posed  in  terms 
of  local  scaling  instead  of  Moni'n-Obukhov  similarity  theory.  However,  the  distinction 
m  the  literature  is  not  always  clear,  since  the  vertical  flux  divergence  is  usually  not 

available  and  the  difference  between  the  flux  at  the  observational  level  and  the  surface 
is  not  known. 
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Appendix 

Elimination  of  records 

The  flux  sampling  errors  and  non-stationarity  for  momentum  and  buoyancy  were 
evaluated  following  Vickers  and  Mahrt  (1997b)  and  Mahrt  (1998).  All  vertical  levels 

°J  dT,f?T™h  one'hour  time  Period  with  flux  sampling  errors  exceeding  prescribed 
thresholds  for  any  of  the  4  levels  for  either  the  momentum  or  the  buoyancy  flux  were 
eliminated  from  further  analysis.  Using  the  notation  in  Vickers  and  Mahrt  (1997b),  the 
thresholds  used  were:  (i)  6  for  extreme  subrecord  flux  outliers  ( event );  (ii)  1  for  the 
relative  systematic  flux  error  (RSE)  associated  with  the  flux  on  scales  larger  than  the 
0-minute  average;  (iii)  0.75  for  the  relative  random  flux  error  due  to  inadequate  sample 
size  (RFE);  and  (iv)  5  for  non-stationarity  of  the  flux  (NR,  defined  in  Mahrt  1998;  his 
Eq.  (1 1)).  The  averaging  time-scale  used  in  the  flux  sampling  error  analysis  was:  (i)  5- 
minutes  for  event  and  RFE,  resulting  in  12  independent  estimates  of  the  flux  in  each  one- 
hour  time  period;  (ii)  5-  and  10-minute  fluxes  for  estimation  of  the  scale-dependence  of 
the  flux  (RSE);  and  (iii)  100-second  fluxes  for  the  non-stationarity  ratio  NR.  All  the  flux 
sampling  criteria  applied  to  the  fluxes  of  momentum  and  buoyancy  for  all  4  vertical 
levels  combined  removed  35%  of  the  one-hour  time  periods.  Excluding  data  based  on 
the  flux  sampling  criteria  reduced  the  scatter  in  the  composited  profiles  of  the  wind  stress 
and  buoyancy  flux.  After  the  flux  sampling  screening,  the  final  flux  dataset  consisted  of 
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415  observations  in  off-shore  flow  and  483  observations  in  on-shore  flow.  No  attempt 
was  made  to  correct  for  potential  flow  distortion,  which  may  contribute  to  scatter  in  the 
relationships. 
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